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1.  STATEMENT  OF  PROBLEM 


This  research  project  covered  a  period  of  approximately  four  years. 
The  program  was  aimed  at  an  understanding  of  the  design  and  fabrication  of 
semiconductor  injection  lasers  for  both  communication  and  high  power 
applications.  Lasers  fabricated  for  communication  systems  are  typically  of 
the  low  power  type  which  emit  power  in  10  to  30  milliwatt  range  while  high 
power  devices  (single  element  types)  have  emission  powers  of  over  50 
milliwatts.  Communication  devices,  fabricated  for  long  wavelength 
(approximately  15,000  Angstroms)  operation  where  optical  fibers  have  both 
low  attenuation  and  dispersion,  use  compounds  of  indium,  gallium,  arsenic, 
and  phosphorus.  On  the  other  hand,  contemporary  high-power  devices  (with 
emission  wavelengths  around  8,000  Angstroms)  are  fabricated  with  aluminum, 
gallium  and  arsenic.  Although  the  research  of  this  project  was  concerned 
with  aluminum  and  gallium  arsenide  compounds,  the  theory  developed  is 
applicable  to  all  semiconductor  lasers. 

The  electromagnetic  resonant  modes  in  a  semiconductor  laser  are 
obtained  from  solving  Maxwell's  equations.  The  ultimate  mode  operation  is 
inextricably  tied  to  the  geometry  of  the  laser  as  well  as  to  its 
interaction  or  coupling  with  the  external  world.  The  particular  mode  that 
operates  is  determined  by  its  interaction  with  the  active  region  of  the 
laser.  Most  system  applications  require  devices  that  emit  power  in  a 
single  and  coherent  beam  of  radiation.  However,  a  majority  of  lasers 
radiate  in  either  multiple  beams  or  a  single  unstable  beam  (drive 
sensitive).  On  the  other  hand,  gas  lasers  can  be  fabricated  to  emit  power 
in  a  single  stable  beam  of  radiation  because  the  modes  in  a  gas  laser  are 


formed  by  the  cavity  mirrors;  the  modal  selection  process  is  thus 
independent  of  drive.  In  an  injection  laser  the  drive  dependent  gain 
region  plays  a  major  role  in  the  formation  of  the  waveguide  geometry. 
Consequently,  we  must  carefully  design  structures  that  can  operate  only  in 
a  single  mode. 

Because  of  design  restrictions,  many  contemporary  lasers  are 
constrained  to  low  powers.  To  increase  the  output  power,  several  exotic 
structures  have  been  proposed  in  the  literature.  But,  these  newly 
conceived  devices  have  pushed  the  fabrication  technology  to  its  limits 
which  adversely  affects  manufacturing  yields. 

Design  sophistication  of  single  laser  elements  is  one  method  of 
achieving  high  overall  emission  powers.  Another  approach  is  to  couple 
several  lasers  together  to  form  a  linear  array.  Since  the  devices  are 
fabricated  on  a  single  wafer,  locking  of  the  elements  to  obtain  a  single 
coherent  source  can  be  accomplished  via  internal  field  coupling.  This 
coupling  is  possible  because  the  waveguides  are  open  types.  Nevertheless, 
proper  design  of  the  single  element  to  achieve  the  "right"  coupling  is 
necessary. 

Thus,  the  problem  of  analyzing  the  individual  laser  element,  the 
central  theme  of  this  research  project,  is  the  single  most  important  aspect 
required  for  the  ultimate  realization  of  the  semiconductor  injection  laser. 


SUMMARY 


INTRODUCTION 

Semiconductor  lasers  are  finding  applications  in  many  areas,  such  as 
optical  recording,  high-speed  printing,  fiber  optic  systems,  communication 
through  free  space,  optical  pumping  of  various  lasers  which  operate  at 
different  wavelength  and  others.  Development  of  high-power  lasers  has  made 
it  possible  to  fabricate  single  element  lasers  that  radiate  over  100 
milliwatts  in  continuous  operation  and  much  higher  powers  in  pulsed 
operation.  Also,  recent  work  in  phased-laser-arrays  has  shown  that  it  is 
possible  to  fabricate  individual  devices  with  relatively  high  packing 
densities.  Thus,  the  overall  power  emission  from  the  single  coherent 
source  is  very  high. 

Semiconductor  lasers  will  begin  to  replace  many  of  the  very  high- 
power  optical  sources  in  the  next  5  to  10  years.  There  are  two  main 
reasons  that  will  propel  the  growth  in  the  high  power  applications.  (1) 
Semiconductor  devices  are  superior  to  "vacuum  tube"  components  in  terms  of 
reliability.  (2)  For  a  given  volume,  the  power  emission  from  a 
semiconductor  laser  compared  to  that  obtained  from  a  gas  laser  is  extremely 
large.  High  powers  are  possible  because  of  the  packing  density  of  excited 
electrons  or  inversion  levels  is  semiconductors. 

The  operation  of  injection  lasers  can  be  categorized  into  several 
areas.  These  include  (1)  power  emission  capacity,  (2)  operational 
wavelength,  (3)  wavelength  or  mode  stability,  and  (4)  spatial  radiation 
stability.  Typically,  many  of  the  above  are  interconnected.  For  example, 
high-power  lasers  tend  to  develop  many  operational  instabilities.  These 


instabilities  occur  because  of  the  strong  nonlinear  relationships  between 
optical  waveguide  formation  in  the  material  and  the  optical  mode  intensity 

Much  of  the  work  accomplished  from  this  grant  concerned  the 
development  and  understanding  of  single-component  high-power  injection 
lasers.  Our  models  of  electromagnetic  phenomena  in  nonplanar  lasers  have 
helped  researchers  in  the  laser  field  to  push  toward  optimum  designs.  At 
this  point  in  time,  many  possible  structures  can  produce  high  powers. 

These  include  the  channel-substrate-planar  (CSP),  the  constricted-double¬ 
heterojunction  large-optical-cavity  (CDH-LOC),  and  others.  In  terms  of 
modeling  the  various  devices,  the  CDH-LOC  laser  is  the  most  difficult.  We 
discuss  some  of  its  properties  below.  Extension  of  these  models  to  other 
nonplanar  lasers  such  as  the  CSP  is  easily  accomplished. 

MODELS  OP  NONPLANAR  SEMICONDUCTOR  LASERS 

I  2 

The  CDH-LOC  laser  provides  '  high  single-mode  power  because  of  the 

nonplanarity  of  the  various  epitaxially  grown  layers.  Because  this  device 

is  nonplanar,  it  is  very  difficult  to  model  many  of  its  subtle 

characteristics.  However,  our  model  of  the  CDH-LOC  optical  modes  is 

confirmed  by  excellent  fits  of  theoretically  generated  far-field  patterns 

to  experimental  ones.  A  major  finding  discovered  by  our  model  was  the 

combined  effect  of  lateral  variations  in  transverse  (optical)  confinement 
3 

factor,  r,  and  gain-induced  index  depressions  on  lateral  mode  control. 

Figure  1  depicts  schematically  the  basic  geometry  of  CDH-LOC  devices 

a  "convex-lens-shaped"  active  layer  atop  a  "concave-lens-shaped"  guide 

layer.  The  bulk  refractive  indices  are  such  that  ni>n^>n^>n^.  We  have 

shown4  that  such  structures  translate  into  passive  W-shaped  lateral 

2 

waveguides.  Unlike  the  waveguiding  structures  usually  analyzed,  the  W- 


guide  has  no  known  solutions  of  the  wave  equation.  Thus  it  was  necessary 
to  employ  numerical  methods  to  solve  Maxwell's  equations  for  arbitrary 
dielectric  profiles,  and  thus  to  solve  for  the  lateral  mode  content  of  the 
W-guides  as  described  below. 


Fig.  1  Cross  section  of  a  constricted-double-heterojunction  laser  with 
a  large-optical-cavity.  Light  is  confined  to  the  regions  1  and 
2  and  to  the  lateral  region  near  y  ■=  0. 


Solutions  to  Maxwell's  equations  were  restricted  to  fields  polarized 
along  the  junction  with  E(x,y)  =  f(x,y)  h(y)  expfjwt  -  yz)  where  f(x,y)  and 
h(y)  describes  the  transverse  and  lateral  fields  respectively.  The  lateral 
modes  are  obtained  from  a  solution  of 

d2h/dy2  =  [y2  -  y2  *  k2  T(y)  ^(yjjh  *  0  (1) 

where  kQ  *  2v/\Q,  with  the  free-space  wavelength  and  a^fy)  is  the 
current  dependent  change  in  the  complex  dielectric  constant  of  the  active 
layer.  The  passive  structures  have  6*^  *  0.  yQ  is  the  one-dimensional 
propagation  constant  that  determines  the  effective  index  of  the  passive 
structure,  n  (n _ *  jy n/kA);  and  y  is  the  two-dimensional  propagation 

CO  fO  U  v 
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constant  to  be  determined.  The  values  of  are  found  from  the  current 
distribution  flowing  into  the  active  layer.  An  estimate  of  the  current 
spreading  was  obtained  by  numerically  solving  an  approximate  circuit 
model. 5,6  We  found  the  current  to  be  virtually  flat  under  the  10  um-wide 
contact  stripe,  and  tapering  to  half  maximum  at  approximately  5  urn  from  the 
contact  stripe  edges.  Using  the  current  distribution  in  the  active  layer, 

7 

we  can  write  the  gain 

g(y)  =  4 5  Jrff (y)/d(y)~  190  (cm'1)  /(2) 

2 

where  J  e{{(kA/ca  )  is  the  effective  current  density,  and  d(y)  is  the  active 

layer  thickness.  The  gain  affects  both  the  imaginary  and  real  parts  of  the 

2 

dielectric  constant  : 

*  2n1«n1(yJ  +  j6r>1  g(y )/kQ  (3) 

where  Sn^y)  is  the  gain-induced  index  change  in  the  active  layer.  (It  is 
proportional  to  the  carrier  injection.)  Since  fin^fy)  is  generally  found  to 
be  negative  we  refer  to  it  as  gain-induced  index  depression.  The  net 
(real)  effective  index  that  characterizes  the  lateral  field  is 

neff (y)  *  neo(y)  +  T(y)  an^y)  (4) 

If  we  assume  that  the  background  absorption  coefficient  aQ  in  the 
unpumped  active  regions  is  190  cm-'1  we  have  for  Sn^  the  following  ad  hoc 
relationship 

=  R  (g(y)  +  190)  cm-1  /kQ  (5) 

with  R  being  the  ratio  of  index  changes  to  normalized  gain  changes,  also 

8  9 

known  as  antiguiding  factor.  '  Finally,  by  using  (2)  and  (5)  we  obtain 

•(6) 


Sn.(y)  =  45  R  J^Ay )/k.d(y) 


The  analysis  was  applied  to  two  different  high-power  CDH-LOC 

structures.  Active  and  guide-layer  thicknesses  were  accurately  measured  in 

the  lateral  direction,  over  the  range  10  vm<y <10  tin,  by  using  5°  angle- 

lapped  stained  cross-sections.  The  A1  content  of  the  various  A1  Ga,  As 

x  1-x 

layers  was  determined  from  photoluminescence  data.  The  layers  varied 

laterally  in  thickness  as  shown  in  Fig.  1  (For  one  geometry  the  guide  layer 

was  asymmetric  due  to  substrate  misorlentation  effects  on  growth. With 

the  structures  fully  characterized,  the  confinement  factor  r(y)  and  the 

effective  index  n  (y)  are  determined  by  solving  for  the  transverse  modes 
60 

at  each  lateral  data  point.  The  y)  plots  for  the  two  structures  are 
shown  in  Figs.  2  (a)  and  (c).  In  the  same  figures  we  plot  the  respective 
net  effective  refractive  index  profiles,  neff(y)  (i.e.f  W-guides).  Both 
the  passive  contribution  n  ( y)  and  the  active  contribution,  r(y)  «n.(y), 

60  X 

are  considered.  While  and  r  were  known,  dn^fy)  had  to  be  deduced  from 

fitting  of  theoretical  far-field  patterns  to  experimental  ones  (see  Figs.  2 

(b)  and  (d) ) .  Thus  we  obtain  for  the  two  structures  at  threshold  dn^(0) 

values  of  -0.028  and  -0.02,  which  agree  with  dn^  values  determined  by  other 
11  12 

workers  '  from  different  experiments.  The  gain  in  the  active  layer  at 
threshold  was  set  at  50  cm  *  to  offset  end  and  internal  cavity  losses. 
Having  generated  the  proper  W-guides,  we  then  solved  the  wave  equation 
numerically  for  the  lateral  modes.  Radiation  patterns  were  determined  from 
the  Fourier  Transform  of  the  near  fields. 

As  far  as  the  discrimination  against  high-order  mode  oscillation,  two 
features  of  Figs.  2  (a)  and  (c)  are  relevant:  1)  r(y)  peaks  at  y  =  0  and 
drops  relatively  fast  across  the  lasing  region  (5-7  vm  wide);  and  2) 
neff*y*  Profi^es  are  W-guides.  Both  features  favor  oscillation  in  the 
fundamental  mode  over  the  first-order  mode.  It  should  be  stressed  that  the 


major  effect  of  the  gain-induced  index  depressions  is  to  change  the  n  ^ 
profile  form  a  mild  W-guide  to  a  strong  one.  This  result  follows  because 
the  contribution  of  the  gain-induced  index  depressions,  6n^( y) ,  to  the  net 
effective  index  (Eq.  (4))  is  weighted  by  the  strongly  varying  r(y) 
function.  Thus,  upon  applying  gain  the  center  of  the  "built-in"  W-guide 
is  depressed  more  than  its  wings.  In  turn  high-order  modes  are  further 
suppressed.  Since  the  r(y)  profile  is  virtually  drive  independent  it 
follows  that  in  CDH-LOC  devices,  we  have  a  stable  waveguide  mode.  By 
contrast,  planar-stripe  lasers  do  not  have  a  built-in  lateral  variation  in 
r,  and  thus  are  vulnerable  to  drive. 

The  far-fields  were  fitted  by  adjusting  only  one  parameter:  R(see 
Eq.  (5).  For  run  DB-208  (Fig.  2(b))  side  lobes  at  10°,  indicative  of 
radiative  leakage,  could  only  be  fit  for  R  =-4.  Run  DB-181  (Fig.  2d)  had  no 
leakage  for  the  fundamental  mode  since  the  central  region  of  the  W-guide 
was  wide.  Yet  a  lateral  beamwidth  of  approximately  6°  could  be  fitted  only 
for  R  values  in  the  -3  to  -4  range.  Once  R  was  determined,  was 
obtained  from  (6),  which  in  turn  allows  the  generation  of  the  actual  W- 
guides . 

The  lateral  modal  content  of  the  two  structures  for  R  values  of 
practical  interest  and  threshold  mode  gain  of  50  cm  *  is  as  follows:  only 
the  fundamental  mode  for  run  DB-208;  and  three  modes  for  run  DB-181. 

Figure  2  (e)  shows  the  threshold  current  density  for  the  first  two  modes  of 
DB-181  lasers  as  a  function  of  the  antiguiding  factor  R.  We  also  plot  the 
index  depression  at  the  center  of  the  waveguide  for  various  and  R 

values.  Note  that  the  fundamental  mode  (i.e.,  Oth  order  mode)  has  a 
relatively  constant  threshold  current  density  for  all  R  values  of  practical 


discontinuity  at  R  =  -2.9.  Above  that  value  the  first  order  mode  is  being 
"pushed  out"  of  the  central  index  guiding  region.  Thus,  the  current 
density  must  be  increased  to  large  values  in  order  to  re-establish  the 
mode.  From  a  practical  stand  point  we  can  consider  the  first-order  mode  to 
be  cut  off  at  R  =  -2.9.  It  is  thus  apparent  that  by  varying  the  gain- 
induced  index  depression  one  can  choose  between  multi-mode  and  fundamental - 
mode  operation.  Changes  in  n ^  can  be  incurred  via  changes  in  device  length, 
facet (s)  reflectivity  and/or  ambient  temperature. 


REFRACT  ive  index  (effective)  refractive  index  (effective) 
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Fig.  2  Optical  modes  of  CDH-LOC  laser  structures  are  calculated  by  the  effective 
index  method,  (a)  and  (c)  are  for  two  different  lasers,  (b)  and  (d) 
show  the  experimental  and  theoretical  far-field  patterns.  "Mode  cutoff" 
conditions  for  DB  181  can  be  estimated  from  (e) . 
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The  research  carried  out  over  this  past  period  involved 
the  exploration  of  the  field  solutions  obtained  by  a  rigorous 
method  developed  by  the  Principal  Investigator.  A  highlight 
of  our  work  shows  that  stripe  geometry  laser  modes  have 
characteristics  that  are  similar  to  leaky  waveguide  modes. 

In  our  previous  work  we  assumed  that  the  dielectric 
profile  in  the  injection  laser  had  a  parabolic  functional 
variation  along  the  lateral  direction(parallel  to  the 
junction  plane) .  Although  the  parabolic  profile  has  been 
extensively  used  by  previous  investigators,  it  has  definite 
realistic  limitations.  In  our  present  work  we  are 
developing  a  model  using  a  dielectric  profile  which  can 
more  accurately  characterize  the  true  one.  For  example, 
at  large  lateral  distances  (away  from  the  current 
injecting  contact)  the  dielectric  constant  must  approach 
a  limiting  value;  in  the  parabolic  profile  the  dielectric 
constant  continues  increasing/decreasing  dramatically 
with  lateral  distances. 
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HKIl'F  OITLINF  Of-  KI.SFAKCH  KINOINCS 

The  work  during  this  past  period  has  been  aimed  at  an 
understanding  of  the  lateral  mode  structure  in  CDH-LOC  and 
stripe  geometry  lasers.  The  CDH-LOC  devices  are  basic  high 
power  lasers  with  large  optical  cavities.  Important  cavity 
dimensions  are:  (1)  thickness  between  heterojunctions  that 
confine  the  light  and  (2)  the  width  of  the  active  region  which 
is  dictated  by  the  injection  current  spread.  The  important 
results  of  the  study  during  this  period  are  that  structures 
can  be  designed  with  varying  cavity  thicknesses  aimed  at 
the  control  of  lateral  modes.  Ordinarily,  lateral  modes  are 
not  affected  by  cavity  thicknesses. 

We  have  had  several  verbal  inquiries  concerning  the  results 
of  our  recent  publication  in  the  Journal  of  Quantum  Electronics . 
This  publication  is  the  first  to  predict  how  lateral  mode 
structures  are  affected  by  thin  optical  waveguides.  Previous 
publications  have  linked  radiation  pattern  sidelobes  to  leaky 
modes  while  we  have  shown  these  lobes  result  from  a  "mode 
mixture",  a  requirement  for  a  proper  boundary  value  solution 
to  the  wave  equation. 


BR IFF  OUTLINE  OF  RESEARCH  FINDINGS 


The  work  in  this  past  period  has  been  directed  toward  the  extension  of 
lateral  mode  analysis  in  constricted  double  heterojunction-large-optical 
(CDH-LOC)  devices.  These  devices  are  fabricated  in  such  a  fashion  that 
the  active  region  as  well  as  the  optical  mode  confining  region  vary  in 
thickness  as  a  function  of  lateral  position.  These  lasers  are  far  superior 
in  performance  when  compared  to  other  injection  laser  devices.  CDH-LOC 
lasers  have  excellent  power  output  and  spectral  purity.  (Device  powers 
around  100  mwatts  radiate  in  a  single  mode.)  The  Principal  Investigator 
is  presently  preparing  the  analysis  of  CDH-LOC  lasers  in  conjunction  with 
the  experimental  program  at  RCA  Laboratories. 

Concurrent  with  the  above  work,  we  are  extending  the  analysis  of  the 
exact  mathematical  solution  of  two  dimensional  modes  in  dielectric  wave¬ 
guides.  It  should  be  mentioned  that  there  are  published  results  for  exact 
mode  solutions  only  for  (1)  optical  fibers,  (2)  dielectric  waveguides  of 
rectangular  cross-section,  and  (3)  slab  waveguides  with  localized  lateral 
gain  distributions.  (The  latter  work  was  previously  reported  by  the  PI.) 
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BRIEF  OUTLINE  OF  RESEARCH  FINDINGS 


During  this  past  six  months  much  effort  has  been  spent  on  the 
development  of  :  (1)  an  understanding  of  leaky  wave  structures  found 

in  contemporary  high-power  injection  lasers,  and  (2)  the  integrity 
of  the  effective  index  method  which  is  widely  used  in  integrated  optics. 

The  work  reported  in  the  first  paper  listed  above  contains 
extensive  information  on  the  analysis  of  leaky  mode  lasers.  This  work 
is  the  first  to  report  and  to  properly  analyze  devices  lasing  in  a 
true  leaky  mode.  It  was  previously  believed  that  leaky  modes  could 
not  reach  threshold  because  of  their  lossy  nature. 

The  work  reported  in  the  second  publication  compares  the  modal 
solutions  of  a  class  of  dielectric  waveguides,  obtained  by  exact  and 
effective  index  methods.  The  waveguides  studied  are  those  which 
closely  approximate  actual  guides  designed  for  injection  lasers.  Our 
conclusion  is  that  the  effective  index  method  is  excellent  for  waveguide 
analysis. 


MODE  CHARACTERISTICS  OF  NON -PLANAR  DOUBLE-HETEROJUNCTION 
AND  LARGE-OPTICAL-CAVITY  LASER  STRUCTURES 


J.  K.  Butler*  and  D.  Botez 
RCA  Laboratories 
Princeton,  NJ  08540 


Abstract 


Mode  behavior  of  non-planar  double-hetero junction  (DH)  and  large- 
optical-  cavity  (LOC)  lasers  is  investigated  using  the  effective  index 
method  to  model  the  lateral  field  distribution.  The  thickness  variations 
of  various  layers  for  the  devices  discussed  are  correlated  with  the  growth 
characteristics  of  liquid-phase  epitaxy  over  topographical  features  (chan¬ 
nels,  mesas)  etched  into  the  substrate.  The  effective  dielectric  profiles 
of  constricted  double-hetero junction  (CDH)-LOC  lasers  show  a  strong  influ¬ 
ence  on  transverse  mode  operation:  the  fundamental  transverse  mode  (i.e., 
in  the  plane  perpendicular  to  the  junction)  may  be  laterally  index-guided 
while  the  first  (high)-order  mode  is  laterally  index-antiguided.  The  ana¬ 
lytical  model  developed  uses  a  smoothly  varying  hyperbolic  cosine  distri¬ 
bution  to  characterize  lateral  index  variations.  The  waveguide  model  is 
applied  to  several  lasers  to  illustrate  conditions  necessary  to  convert 
leaky  modes  to  trapped  ones  via  the  active-region  gain  distribution. 
Theoretical  radiation  patterns  are  calculated  using  model  parameters,  and 
matched  to  an  experimental  far-field  pattern. 

This  work  has  been  partially  supported  by  the  U.  S.  Army  Research  Office 
under  Grant  Ho.  DAAG29-80-K-004  and  by  NASA,  Langley  Research  Center, 
Hampton,  VA  under  Contract  No.  NAS1-15440. 
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Comparison  of  Numerical  and  Effective  Index  Methods 
for  a  Class  of  Dielectric  Waveguides* 

A.  Linz  and  J.  Butler 
Electrical  Engineering  Department 
Southern  Methodist  University 
Dallas,  Texas  75275 


Abstract 

A  numerical  method  and  the  effective-index  method  are  applied 
to  a  three-layer  dielectric  waveguide  with  Eckart  type  dielectric 
constant  variation  in  the  active  layer.  The  results  of  the  2  methods 
are  compared  in  terms  of  the  propagation  constant  y. 
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BRIEF  OUTLINE  OF  RESEARCH  FINDINGS 


In  the  work  over  the  last  period  we  have  continued  the  investigation 
of  approximate  methods  used  for  the  analysis  of  mode  characteristics  of 
injection  laser  devices.  We  have  reported  that  contemporary  Injection 
lasers  with  wide  active  regions  can  be  accurately  modelled  with  the 
effective  index  method.  This  is  a  very  important  result  since  many  new 
lasers  have  non-planar  layers;  the  effective  index  method  affords  a 
convenient  mode  analysis  procedure.  It  should  be  noted  however,  that  there 
are  no  exact  methods  for  analyzing  modes  of  non-planar  laser  structures. 

It  is  assumed  that  since  approximate  and  exact  calculations  of  mode  character 
on  planar  structures  produce  almost  Identical  results,  that  the  approximate 
solutions  are  satisfactory  for  non-planar  devices. 

In  our  last  publication  we  have  theoretically  predicted  transverse 
mode  operation  of  a  multimode  laser.  Mode  operation  was  predicted  by 
calculation  of  the  coupling  between  the  optical  field  and  injected  charge 
density.  We  experimentally  investigated  mode  switching  by  heating  a  diode 
laser;  at  low  temperatures  the  device  operated  In  the  first  high  order 
mode  while  at  high  temperatures ,  fundamental  mode  operation  was  observed. 

The  prediction  of  switching  was  tied  to  the  fact  that  the  optical  field/ 
charge  distribution  coupling  are  temperature  sensitive. 
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Dallas,  Texas  7527S 

Abstract 

Approximate  methods  are  used  to  obtain  the  modal  properties 
of  stripe-contact  semiconductor  injection  lasers  using  a  planar 
three-layer  waveguide  model.  The  central  active  layer  has  a 
dielectric  constant  that  varies  smoothly  along  the  direction 
parallel  to  the  heterojunction  boundaries.  The  complex  dielectric 
constant  under  the  stripe  contact  is  dependent  on  the  gain  and 
approaches  a  constant  value  at  large  lateral  distances.  The  two 
methods  are  compared  in  terms  of  their  modal  propagation  constants. 
An  application  of  the  effective  index  method  facilitates  a  physical 
understanding  of  dielectric  waveguide  modes  as  well  as  providing  an 
efficient  calculation  procedure. 
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ABSTRACT 

Transverse  mode  selection  is  characterized  for  GaAs/AlGaAs  double 
heterojunction  lasers  from  optical  field  and  electron/hole  interaction. 

The  electron/hole  distribution  determined  from  a  solution  of  the  ambipolar 
diffusion  equation  provides  the  necessary  information  about  gain/mode 
coupling  to  predict  the  current  at  threshold.  Lasing  power  out  versus 
current  solutions  provide  information  about  internal  differential  quantum 
efficiency.  Theory  is  matched  to  experiment  for  a  multimode  laser  with  one 
heterojunction  having  a  very  small  index  step.  It  is  found  that  the  laser's 
characteristics  over  a  temperature  and  current  range  are  predicted  by  adjusting 
the  active  layer  refractive  index  as  determined  from  far  field  measurements. 
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BRIi  1  OUTLINE  OF  RESEARCH  FINDINGS 


The  work  over  this  last  period  has  been  concerned  with  both 
the  theoretical  and  experimental  characterization  of  high-power 
semiconductor  laser  devices.  These  lasers  will  have  extensive 
applications  in  military  communication  and  data  processing  sys¬ 
tems.  The  particular  devices  modelled  are  the  constricted- 
double-hetero junction  large-optical-cavity  (CDH-LOC)  lasers  which 
are  fabricated  at  RCA  Laboratories.  These  lasers  have  the 
highest  output  power  (approximately  50  mW  in  continuous 
operation)  of  any  contemporary  injection  laser.  Furthermore, 
their  operational  range  for  single  mode  lasing  is  excellent.  Our 
laboratory  has  developed  the  most  comprehensive  theoretical  model 
of  the  modal  characteristics  of  CDH-LOC  devices. 

Typically,  reasonable  theoretical/experimental  performance 
comparisons  are  made  via  the  far  field  radiation  patterns.  To 
augment  our  theory  we  have  developed  computer  controlled 
apparatus  for  taking  far-field  data  in  two  dimensions.  For 
example,  the  attached  figure  shows  the  experimental  pattern  of  a 
typical  CDH-LOC  laser  where  8  and  <}>  are  the  two  transverse 
angles  in  the  far  field.  (The  main  lobe  peaks  at  0  =  $  =  0.) 

We  are  now  developing  software  to  transmit  the  experimental 
data  to  a  large  computer  for  the  purpose  of  theoretical/experi¬ 
mental  comparisons.  The  addition  of  these  laboratory  experiments 
to  our  program  will  expedite  our  understanding  of  high-power 
injection  lasers  and  thus  improve  our  overall  theoretical 
modelling  capabilities. 


Modal  Solutions  of  Active  Dielectric  Waveguides 
by  Approximate  Methods 


A.  Linz  and  J.K.  Butler 
Electrical  Engineering  Department 
Southern  Methodist  University 
Dallas,  Texas  75275 


Abstract 


Approximate  methods  are  used  to  obtain  the  modal  properties 
of  stripe-contact  semiconductor  injection  lasers  using  a  planar 
three-layer  waveguide  model.  The  central  active  layer  has  a 
dielectric  constant  that  varies  smoothly  along  the  direction 
parallel  to  the  heteroj unction  boundaries.  The  complex  dielectric 
constant  under  the  stripe  contact  is  dependent  on  the  gain  and 
approaches  a  constant  value  at  large  lateral  distances.  The  two 
methods  are  compared  in  terms  of  their  modal  propagation  constants. 
An  application  of  the  effective  index  method  facilitates  a  physical 
understanding  of  dielectric  waveguide  modes  as  well  as  providing  an 
efficient  calculation  procedure. 


PREDICTION  OF  TRANSVERSE  MODE  SELECTION  IN  DOUBLE  HETEROJUNCTION 
LASERS  BY  AN  AMBIPOLAR  EXCESS  DIFFUSION  SOLUTION 

by 

J.  B.  Delaney 

TRW  Optoelectronics  Division 
R.  R.  Shurtz,  II 

Night  Vision  and  Electro-optics  Laboratory 
J.  K.  Butler 

Southern  Methodist  University 


Abstract 

Transverse-mode  selection  is  characterized  for  GaAs/AlGaAs 
double  heteroj unction  lasers  from  optical  field  and  electron/hole 
interaction.  The  electron/hole  distribution  determined  from  a 
solution  of  the  ambipolar  diffusion  equation  provides  the 
necessary  information  about  gain/mode  coupling  to  predict  the 
current  at  threshold.  Lasing  power  out  versus  current  solutions 
provide  information  about  internal  differential  quantum 
efficiency.  Theory  is  matched  to  experiment  for  a  multimode 
laser  with  one  heteroj unction  having  a  very  small  index  step.  It 
is  found  that  the  laser's  characteristics  over  a  temperature  and 
current  range  are  predicted  by  adjusting  the  active-layer 
refractive  index  as  detemined  from  far-field  measurements. 
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Field  Solutions  for  the  Lateral  Modes  of  Stripe 
Geometry  Injection  Lasers 

JEROME  K.  BUTLER,  member,  ieee,  and  JOSEPH  B.  DELANEY 


Abstract- Thi*  paper  compares  die  lateral  mode  patterns  of  stripe 
geometry  lasers  produced  by  two  theories:  1)  an  approximate  effective 
dielectric  method  and  2)  an  exact  method.  The  exact  method  for  cal¬ 
culating  modes  uses  an  expansion  of  Hermite-Qaussian  functions.  For 
small  normalized  frequencies  and  a  depression  of  the  refractive  index 
under  the  stripe,  the  fundamental  lateral  mode  shows  structure  usually 
associated  with  leaky  waveguide  modes. 

A  COMMON  description  of  the  lateral  modes  of  a  stripe 
geometry  laser  uses  the  Hermite-Gaussian  (HG)  func¬ 
tion  [l]-[8]  while  the  transverse  Held  dependence  is  domi¬ 
nated  by  the  large  index  steps  of  the  heterojunctions  which 
sandwich  the  active  layer.  So  far  only  two  methods  have  been 
proposed  to  tie  together  the  lateral  and  transverse  field  de¬ 
pendencies  under  the  stripe.  The  first  is  the  effective  index 
method  [7]  and  the  second  is  an  exact  method  [8] .  The  ef¬ 
fective  index  method  results  in  a  single  term  lateral  field 
description  with  appropriate  “effective”  parameters  derived  to 
account  for  the  influence  of  the  cladding  layers  on  the  lateral 
mode  shape.  The  advantage  of  this  representation  is  that  a 
single  term  expression  can  easily  be  matched  to  such  experi¬ 
mental  observables  as  near-  /far-fleld  shapes  and  virtual  waist 
positions.  The  exact  description  of  the  lateral  mode  takes 
into  account  the  vertical  geometry  by  matching  the  usual  field 
quantities  at  the  heterojunctions.  The  result  is  a  lateral  mode 
description  by  a  linear  combination  of  HG  functions  with  ap- 
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Fig.  1.  Cross  section  of  a  stripe  geometry  laser. 


propriate  weighting  coefficients.  This  paper  reports  on  a  com¬ 
parison  of  the  two  methods.  A  major  difference  between  the 
two  is  the  divergence  of  the  modal  shape  when  the  lateral  in¬ 
dex  is  depressed  under  the  stripe  and  the  normalized  frequency 
for  the  vertical  geometry  (see  Fig.  1),  v  *  (d/2)  k0  (nj  -  nl)lf2 , 
is  small.  In  general,  the  coefficients  of  the  higher  order  HG 
functions  rise  in  magnitude  as  more  light  escapes  into  the 
transverse  cladding.  For  a  negative  depression  in  the  lateral 
index,  the  higher  order  terms  appear  as  shoulders  on  the  near 
and  far  fields  of  the  fundamental  mode  and  the  modal  shape  is 
deformed  from  the  pure  Gausrian.  We  begin  with  a  brief  re¬ 
capitulation  of  the  wave  equation  solutions  of  each  method. 

The  complex  dielectric  constant  can  be  written  as  [7]  (see 
Fig.  1) 

[x>(x)  =  x0-  fro***1.  Ij'xf 

k(*,.vW  d  0) 

“«c.  \y\>2 

where  x0  is  the  complex  dielectric  constant  at  x  »  0,  k®  is  the 
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free  space  wavenumber,  and  kc  are  the  dielectric  constants 
of  the  cladding  layers.  The  parameter  “a”  in  (1)  is  defined  as 


-  (2*  +  /10-4/*o),/1 

k0x0 


a 

R-  —  cm. 


(2) 

(3) 


x0  is  defined  as  the  point  at  which  gain  has  fallen  to  zero,  i.e., 
6 g  -  g0.  6n  is  the  refractive  index  change  at  x  ■  x0.  In  the  ef¬ 
fective  index  method,  the  wave  equation 

(Jfcjic  ( x,y )  +  7s]  ♦  * 0 

with 


♦  «  4/(x)  <p  O')  e'7* 


dy 2  ' 


»>»<! 

\y\>  | 

(4) 

(5) 

(6) 


is  multiplied  by  4'*  and  integrated  over  ye(~°°,  “>)  to  yield  an 
effective  solution 

*2x0r+^x«(i  -  r>-  ^r+p*(i  -  d+72 
-*5«r,'a(2/+i) 

«.rf  “  r**« 

where  T  is  the  confinement  factor 

prf/s 


^-<t/2 


dO,)0*O)dy 


(7) 

(8) 
(9) 

(10) 


when  daO)  is  normalized  to  unity,  and 
(cos  qy 


t(y)mA 


ljl<f 


|COffl|eP[W/s)-b'U  „|>|. 


01) 


The  quantities  p  and  q  are  the  eigenvalues  of  the  slab  wave¬ 
guide.  In  17] ,  the  equation  corresponding  to  our  (7)  did  not 
include  the  eigenvalue  p.  The  integration  over  the  y  direction 
dicutes  an  effective  represenution  of  p*(l  -  T)  -  q*  r  u  in¬ 
cluded  in  (7). 

The  exact  method  can  be  summarized  [8] 

♦  ■  ^  At  cos  (qfy)  Hi  (fl1/s**x)  e’,/***»*‘ 


\y\<\ 

•J  ^0f)«osX*«*p|^"3^  (X*  "  7a)1/,|^X 

lyl>|  (12) 


ACTIVE  REGION  THICKNESS  d  (micron,) 

NORMALIZED  FREQUENCY  • 

FIs.  2.  Far  Bald  half-power  hill  width  #|  versus  nonnalised  frequency  v. 

where  /  ■  0, 2, 4,  ■  •  •  for  even  modes.  The  vertical  eigenvalue 
qi  satisfies 

7*  ««?-*&* o  +  k|««(2I+l).  (13) 

Parameters  used  in  sample  calculations  are  ngmnem  3.4, 
nb( 0)  ■  3.6  X  *  05  pm,  G  ■  50  cm-1,  and  o,  ■  ae  *  20  cm* *. 
Comparison  of  the  two  models  begins  with  the  gain.  For  the 
effective  index  method  (7)  can  be  used.  The  term  ko(jc0r  + 
x«(l  -  T))  is  the  effective  dielectric  constant,  while pJ(l  -  T)- 
represents  an  effective  vertical  eigenvalue.  Vrtth  y  * 
-  C/2  +  tfi,  the  gain  in  the  active  region  must  include  the  mirror 
loss,  losses  in  the  cladding  layers,  and  losses  in  the  active  layer 
represented  by  the  rigit-hand  side  of  (7).  For  instance,  with 
R  <  0,  the  gain  region  requires  additional  pumping  to  account 
for  lateral  refraction  losses.  Since  these  are  the  same  losses  ac¬ 
counted  for  by  the  exact  technique,  both  methods  calculate 
virtually  the  same  peak  gain  values  for  all  values  of  R 
investigated. 

For  R  >  0,  the  near-  and  far-field  structures  coincide  for 
both  methods.  A  peculiarity  of  the  patterns  is  the  behavior  as 
the  normalized  frequency  v  gets  small  (d  -*■  0).  The  far-field 
half-power  full  width  0j  is  plotted  in  Fig.  2  versus  active  region 
width  d.  As  the  active  region  shrinks,  peak  gain  values  increase 
to  satisfy  additional  losses.  The  parameter  ua,”  as  computed 
from  (2),  requires  that  the  near  field  is  pulled  in  and  the  far 
field  spreads  as  the  result  of  high  cavity  gain  values. 

The  two  methods  diverge  in  modal  shape  for  R  <  0  and 
small  normalized  frequencies.  For  the  refractive  indexes  of 
this  work  this  corresponds  to  d<  02  pm.  As  power  leakage 
to  the  cladding  increases,  the  magnitude  of  the  expansion  co¬ 
efficient  At  increases  as  illustrated  in  Fig.  3.  With  R  <  0, 
|«tl>|«r|,  where  amar  +  iat,  and  shoulders  develop  on  the 
fundamental  mode  patterns.  An  example  of  a  far-field  pattern 
evolution  with  d  ■  0.1  pm  is  shown  in  Fig.  4  where  the  dashed 
curve  is  die  fundamental  mode  pattern  for  the  effective 
method  and  the  solid  curve  the  fundamental  mode  pattern  for 
the  exact  method.  The  negative  value  Sn,  computed  from  (3), 
is  a  measure  of  index  change  and  antiguiding  over  the  gain  re¬ 
gion.  It  can  be  seen  that,  fot  a  value  of  8n*-OJ00185  in 
Fig.  4(b),  the  shoulders  on  the  far-field  pattern  are  well  de¬ 
veloped.  Since  the  Fourier  transform  of  an  HG  function  is  an 
HG  function,  the  near-field  modal  structure  is  theoretically 
similar  to  the  shapes  of  Fig.  4.  From  Fig.  4(a)  it  would  appear 
that  there  is  no  difference  in  the  two  models  for  R  ■  -10'* 
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NORMALIZED  FREQUENCY  v 

FI*.  3.  Magnitude  of  the  expansion  coefficient  Aj  remit  normalized 
frequency  v.  At  more  light  spreads  vertically  into  the  cladding  layers, 
e  larger  percentage  of  higher  order  Hermite-Geunias  functions  are 
needed  to  match  appropriate  field*  at  the  heterojunctions. 


(«)  <b)  (c) 

Fig.  4.  Far-Held  lateral  intensity  of  the  fundamental  mode.  The  dashed 
lines  ate  for  the  effective  method.  For  all  curvet  d  ■  0.1  Mm, 
numnem  3.4, n* “  3.4. 


and  d  ■  0.1  Jim.  In  fact,  since  the  modal  shapes  are  so  similar 
a  curve  for  this  case  R~- 10"*  is  included  in  Fig.  2.  A  final  re¬ 
mark  about  Fig.  4  is  that  the  position  of  the  shoulders  is  pre¬ 
dominately  determined  by  the  gain  width  2x0  and  not  R.  An 
index  depression  under  die  stripe,  represented  by  a  negative 
value  of  R,  determines  the  intensity  of  the  shoulders  but  not 
the  lateral  position. 

In  conclusion,  we  have  briefly  compared  two  solutions  for  a 
parabolic  lateral  dielectric  constant  in  the  active  layer  of  a 
stripe  geometry  laser.  Both  solutions  assume  that  the  vertical 
dependence  is  due  to  die  well  defined  index  steps  at  the 
heterojunctions.  The  effective  and  exact  methods  predict  the 
same  peak  gain  values  since  both  satisfy  the  same  loss  mecha¬ 
nisms.  The  two  methods  diverge  in  modal  shape  for  the  case 
of  a  depression  in  the  refractive  index  in  the  gain  region  below 
the  stripe  and  small  normalized  frequencies.  Matching  the 
fields  at  the  vertical  heterojunctions  results  in  a  fundamental 
modal  shape  that  has  shoulders  on  both  the  near-  and  far-field 
patterns  in  the  case  of  lateral  antiindexing  guiding. 
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Of  primary  Importance  frr  many  applications  h  the 
amount  of  ew  power  that  can  ba  delivered  rcHably 
bom  a  dtoda  laser  Into  a  alngle  mode  (spatte'  -.id  In 
frequency).  One  approach  toward  *te  ach'ivpment 
of  Mgfvpower  mode-stabilized  diodes  It  the  ore  of 
large -opttcal-cavlty  (IOC)  struetrass.'  *  frwt  la, 
» true  lure*  In  which  the  optical  mode  acquires  gain 
from  a  IWn  (0.05-0.2-pm)  active  layer,  while  a  siz¬ 
able  part  of  the  mode  energy  propagates  In  tlilck 
(0.5-I.S-jim)  guide  layedt)  adfeceot  to  the  active 
layer.  We  have  reported  previously  on  hl^vpower 
leaky-cavtty  COH-LOC  lasers.1  Such  de vices 
provide  single-mode  cw  operation  up  to  40  mW  from 
one  facet  In  beams  of  narrow  transverse  full  width 
(8±  "  25-30*)  and  In  large  lasing  spots  (1.5  X  • 
pm).  Here  we  report  on  a  new  type  of  COH  device: 
positive- Index  COH-LOC.  and  we  present  a  com¬ 
prehensive  treatment  of  lasing  mode  confinement 
and  selection  In  COH-LOC  structures.  Positive- 
Index  COH-LOC  devices  have  a  relatively  thick 
(0  20-0.35-pm)  convex-lens-shaped  active  layer 
grown  above  a  concave-lens-shaped  guide  layer. 
By  comparison  with  leaky-cavtty  devices,  (t)  no 
radiation  side  lobes  are  observed  In  the  lateral  far- 
fleld  pattern  (#|  ”  8*):  (2)  the  transverse  far  field  Is 
stilt  relatively  narrow  8L  ~  30s);  (3)  the  threshold 
currents  are  lower  (80-70  mA  vs  80- 130  mA)  for 
the  same  device  length;  and  (4)  single-mode  oper¬ 
ation  Is  achieved  to  only  12  mW/facet  (cw)  and  20 
mW/lacet  (pulsed).  Also  the  threshold  curents  are 
found  to  have  record  Mgh-temperature  coefficients 
for  LOC-type  structures;  7*  a  135*0  In  both 
pul  .ed  and  cw  operation,  while  the  difference  In  Al 
concentration  between  the  active  and  guide  layers 
la  only  15%. 

By  using  the  effective-index  method  It  Is  shown 
that  In  LOC  structures  the  combination  ol  a  con¬ 
vex-lens-shaped  active  layer  and  a  concave-lens¬ 
shaped  guide  layer  provides  a  W-sheped  lateral 
waveguide.  Depending  on  the  COH-LOC  structrae 
geometry,  these  W-type  guides  support  relatively 
large  (5-7-pm)  fundamental  optical  modes  that  can 
be  (1)  totally  leaky  (Is  ,  Improper  mode);  (2)  partially 
guided  and  partially  leaky  (Is.,  gulded-leaky  mode); 
or  (3)  totally  guided.  Aside  from  bringing  flexibility 
In  the  COH-LOC  single-mode  laser  design,  W-type 
lateral  guides,  under  certain  conditions,  strongly 
suppress  oscillation  of  both  htgtvorder  lateral  and 
transverse  modes,  much  more  so  than  ridge 
guktes,~*'*7  (local  Increase  h  Index)  or  leaky  guides' 
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(local  decrease  In  Index).  The  traoretlcol  freafrnent 
Is  backed  by  experimental  results  from  COH-LOC 
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Abstract- Mode  behavior  of  nonplanar  doubto-heteroj unction  (DH) 
and  laiRe-optical-cavity  (LOC)  lasera  i*  inveatigated  uainf  the  effective 
index  mediod  to  model  the  lateral  field  distribution.  The  thickness 
variations  of  various  layer*  for  the  devioea  d jammed  are  correlated  with 
the  growth  characteristic*  of  liquid-phase  epitaxy  over  topographical 
features  (channels,  mesas)  etched  into  the  substrate.  The  effective  di¬ 
electric  profiles  of  constricted  double -heterojunction  (CDH>LOC  lasers 
show  a  strong  influence  on  transverse  mode  operation:  the  fundamental 
transverse  mode  (i-*-,  in  the  plane  perpendicular  to  the  junction)  may 
be  laterally  index-guided,  while  the  Tint  (hjgh>Ofder  mode  is  laterally 
index -antiguided.  The  analytical  model  developed  uses  a  smoothly 
varying  hyperbolic  cosine  distribution  to  characterise  lateral  index 
variations.  The  waveguide  model  is  applied  to  several  lasers  to  illus¬ 
trate  conditions  necessary  to  convert  leaky  modes  to  trapped  ones  via 
the  active-region  gain  distribution.  Theoretical  radiation  patterns  ate 
calculated  using  model  parameters,  and  matched  to  an  experimental 
far-field  pattern. 


I.  Introduction 

ABRICATION  of  semiconductor  injection  laser  diodes 
over  the  past  several  years  has  progressed  to  a  fairly 
sophisticated  level  in  the  sense  that  many  contemporary  de¬ 
vices  are  grown  on  substrate  material  that  has  been  processed 
by  the  incorporation  of  various  grooves  and  undulations  at  the 
substrate  surfaces  [1]-[1S].  After  the  substrate  has  been 
properly  processed,  various  layers  of  GaAs  and  AlGaAs  com¬ 
pounds  are  grown  over  the  grooves.  The  effects  of  grooves  are 
to  produce  active  layers  as  well  as  optical  confining  layers  of 
relatively  small  lateral  thickness  variations,  which  tend  to 
confine  the  optical  radiation  in  the  lateral  direction  defined 
by  the  grown  layers.  In  this  paper  we  discuss  the  performance 
of  devices  with  grown  layers  as  illustrated  in  Figs.  1  and  2. 
Specifically,  we  will  discuss  the  character  of  the  effective 
lateral  refractive  index  and  the  resulting  modes  in  several  types 
of  laser  devices.  Much  of  the  previous  work  on  laser  structures 
has  centered  around  double-heterojunction  devices,  where  the 
active  region  is  sandwiched  between  two  optical  confining 
layers,  as  illustrated  in  Fig.  1(a). .  We  will  discuss  mainly  the 
mode  character  of  the  Helds  in  structures  illustrated  in  Fig. 
1(b)  and  in  high-power  LOC  (large-optical-cavity)  devices  [11], 
where  the  majority  of  the  optical  field  intensity  is  confined  to 
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REGION  REGION 


Fig.  1.  Active  layer  grown  of  nonplanar  double  heterojunction  laaers 
(a)  convex  and  (b)  concave  waveguides. 


Fig.  2.  Nonplanar  large  optical  cavity  (LOC)  User  with  a  GaAs  active 
region  2.  The  optical  cavity  includes  both  Uyers  2  and  3. 

a  region  adjacent  to  the  active  layer  supplying  optical  power 
to  the  lasing  mode.  The  LOC  devices  [10],  [11]  have  been 
reported  to  exhibit  stable  lateral  mode  operation  and  high  out¬ 
put  powers.  Basically,  there  are  two  types  of  LOC  structures 
reported  recently  that  exhibit  stable  lateral  mode  operation. 
These  devices  are  derivatives  of  that  illustrated  in  Fig.  2.  In 
both  cam  the  thickness  of  the  active  layer  decreases  as  a  func¬ 
tion  of  distance  with  the  maximum  value  occurring  in  the 
region  where  die  optical  Held  peaks.  The  adjacent  layer  which 
forms  the  large  optical  cavity,  however,  differs  in  the  two 
structures;  one  [10]  has  the  layer  thickness  decreasing  with 
lateral  distance,  while  the  other  [11],  [IS]  has  the  layer  in¬ 
creasing  with  lateral  distance.  Thus,  the  former  structure  has 
layers  that  act  as  convex  lenses  confining  the  optical  radiation, 
while  the  latter  has  one  layer  causing  convexity  (i.e.,  guiding) 
but  the  other  causing  concavity  (i.e.,  antiguiding). 

After  an  initial  discussion  of  the  waveguide  theory  that  per¬ 
tains  to  nonplanar  lasers,  we  will  discuss  the  lateral  mode 
character  of  lasers  u  they  are  influenced  by  tensing  effects. 
We  address  how  various  layer  thicknesses  of  LOC  devices 
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affect  transverse  and  lateral  mode  operation.  For  example, 
there  exist  combinations  of  thickness  values  dt  and  dJt  which 
determine  whether  the  laser  operates  in  the  fundamental 
transverse  mode  or  the  first  high-order  one.  We  will  present  a 
series  of  design  curves  which  facilitate  the  designs  of  funda¬ 
mental  transverse  mode  lasers.  We  also  find  that  there  exists 
a  peculiarity  that  when  the  laser  operates  in  the  fundamental 
transverse  mode,  the  mode  sees  a  convex  lateral  index  varia¬ 
tion  which,  of  course,  tends  to  confine  the  mode.  On  the 
other  hand,  when  the  laser  operates  in  the  first-order  trans¬ 
verse  mode,  that  mode  sees  a  concave  lateral  refractive  index 
variation,  which  tends  to  defocus  the  mode. 

To  correlate  effective  dielectric  variations  with  practical 
devices  we  present  typical  layer  thickness  variations  in  mode- 
stabilized  DH  and  LOC  structures  grown  on  nonplanar 
substrates.  Our  particular  emphasis  will  be  on  constricted- 
double-heterojunction  (CDH)  [6] ,  [7] ,  [13]  and  constricted- 
double-heterojunction  large-optical  cavity  (CDH-LOC)  [11], 
[IS]  lasers. 

The  lateral  mode  behavior  is  modeled  using  a  hyperbolic 
cosine  variation  of  the  dielectric  constant.  This  analytical 
model  can  be  applied  by-and-large  to  a  broad  number  of  non¬ 
planar  guides.  Another  useful  aspect  is  the  applicability  of  the 
model  to  leaky  modes  which  might  exist  in  index-antiguided 
structures.  The  incorporation  of  lateral  gain  in  the  model 
allows  us  to  predict  when  leaky  modes  become  trapped. 
(True  leaky  modes  have  lateral  fields  that  increase  without 
limit,  while  trapped  modes  have  fields  that  decay  to  zero  at 
infinite  lateral  distances.)  After  developing  fundamental 
characteristics  of  the  lateral  modes  found  from  the  analytical 
field  solutions,  we  apply  our  model  to  specific  CDH  and 
CDH-LOC  lasers.  The  results  indicate  that  most  practical 
devices  lase  in  trapped  waveguide  modes  whether  the  mode  is 
index-guided  or  index-antiguided.  However,  we  address  the 
possibility  that  some  CDH  lasers  [6]  might  oscillate  in  a  leaky 
mode.  The  infinite  field  strength  that  characterizes  leaky 
inodes  never  develops  because  of  the  scattering  of  energy  into 
waveguide  radiation  modes  at  the  facets. 

Finally,  the  radiation  patterns  are  given  for  the  fundamental 
lateral  mode  found  from  analytical  field  solutions.  The  pat¬ 
tern  of  a  CDH-LOC  device  is  measured  and  compared  with  the 
theoretical  pattern  using  the  device  parameters  pertinent  to 
the  hyperbolic  cosine  variation  of  the  dielectric  constant. 

II.  Waveguide  Theory 

The  analysis  of  waveguide  modes  in  contemporary  laser 
structures  is  typically  effected  in  an  approximate  manner. 
Early  approaches  of  analyzing  two-dimensional  waveguides 
other  than  those  encountered  in  optical  fibers  are  to  treat  the 
two  perpendicular  transverse  directions  independently  [16]. 
More  sophisticated  approximate  methods  have  been  refined  to 
what  is  generally  now  recognized  as  the  effective  index 
method  [17],  [19].  In  laser  structures  we  define  two  direc¬ 
tions:  x,  the  direction  perpendicular  to  the  grown  epitaxial 
layers;  and  y,  the  direction  along  the  various  grown  layers  and 
in  the  plane  of  the  mirror  facet.  The  effective  index  method 
yields  Held  solutions  which  accurately  approximate  the  actual 
fields  when  the  transverse  dielectric  variations  are  large  com¬ 


pared  to  the  lateral  ones.  In  stripe  geometry  devices,  the 
lateral  dielectric  dependencies  are  due  to  the  lateral  gain 
variations  which  affect  the  imaginary  part  of  the  dielectric 
constant,  while  temperature  and  carrier  injection  [20]  affect 
the  real  part  of  the  dielectric  constant  (local  temperature  rises 
increase  the  refractive  index  while  carriers  reduce  it). 

The  dielectric  constant  in  the  grown  layers  is  typically  in¬ 
variant  with  respect  to  position.  However,  in  active  devices, 
the  dielectric  constant  of  the  active  layer  (region  2  in  Fig.  1) 
must  reflect  the  gain  distribution.  Then  the  active-layer  dielec¬ 
tric  constant  k2  can  be  written  as 

*7(y)-n\+ig(y)ntlk0  (la) 

=  Koj+Kj„O0  (lb) 

where  n2  is  bulk  refractive  index  of  the  active  layer,  g(y)  is 
the  lateral  gain  distribution,  k0  is  the  free-space  propagation 
constant  ( ka  =  2»A).  end  k02  =  rt\.  It  should  be  noted  that 
when  index  changes  such  as  due  to  temperature  occur,  this 
small  variation  will  be  included  in  k2v(>)  as  a  real  part.  To 
make  bur  analysis  more  general  in  nature,  we  write  the  dielec¬ 
tric  constant  of  the  ith  layer  as 

K/O0  =  *ot  +  k«,/O0,  i  =  1, 2,  3, 4  (2) 

where  nul  contains  only  “lateral"  variations. 

Solutions  to  Maxwell’s  equations  will  be  restricted  to  fields 
polarized  along  the  junction  plane  and  will  be  assumed 
^{x,y)txp{iu3t  -  yz),  where  y  =  a/2  + ip  is  the 
complex  propagation  constant.  The  wave  equation  is 

7?  +  [*ok(x,  y)  +  7*  ]  \l/  =  0  (3) 

where  x(x,y)  defines  the  complex  dielectric  constant  of  all 
space  and  7 ?  is  the  two-dimensjonal  Laplacian  operator.  In 
the  spirit  of  the  effective  index  method,  we  write 

Hx,y)  =  f(x,y)g(y)  (4) 

where  f(x,  y)  describes  the  transverse  fields  and  its  y  depen¬ 
dence  appears  there  because  of  the  layer  thickness  variations. 
The  functional  dependence  of  /  is  determined  from  the  df s 
and  Kol't.  Substituting  (4)  into  (3)  gives 

IfroKOf.^  +  T1]/*'0  (5> 

where  we  have  neglected  derivatives  of  /  with  respect  to  y 
(typically  thickness  variations  iare  relatively  slow  in  the  y 
directions  compared  to  dielectric  variations  in  the  x  direction). 
Equation  (5)  is  now  multiplied  by  f  *  and  integrated  on  x  in 
the  range  (-•»,  °°).  This  yields 

0 +  y'g  +  -  r,o)h?oo)J  g = o  (6) 

where 

*?0)  ■  “  J  (6») 

and  the  overlap  parameter  rV(y)  satisfies 

r,*  f  l/(x,y)l3  dx.  (6b) 

•'fth  layer 


3.525 


The  transverse  functions  f(x,y )  are  normalized  to  unity  along 
x  so  that  I/P/  =  1.  In  the  effective  index  approach  a  propaga¬ 
tion  constant  70,  which  appears  in  a  two-dimensional  wave¬ 
guide  model,  satisfies 

H*h}-k\K0i.  (7) 


Substituting  (7)  into  (6)  gives  the  differential  equation  describ¬ 
ing  the  lateral  fields 


*2± 

dy3 


+ 


(8) 


which  is  equivalent  to  the  result  derived  by  a  more  rigorous 
method  (19].  It  should  be  noted  that  for  layered  waveguides 
with  no  ohmic  losses,  the  effective  propagation  constant 
y0  ■  i0o  is  imaginary  for  all  proper  modes.  At  this  point  we 
should  note  that  there  are  mainly  two  types  of  contemporary 
lasers  considered  here:  1)  double-heterojunction  devices  which 
have  dj  =  0  and  dt  —  0.1 -0.3  pm,  and  2)  LOC  structures  with 
</]  as  above  and  with  d,  ~  1-2  pm.  Ordinarily,  the  former 
structure  supports  only  a  fundamental  transverse  “proper” 
(i.e.,  trapped)  mode,  while  the  latter  can  support  high-order 
proper  modes.  In  LOC  devices  the  relative  values  of  d7  and  d3 
should  be  chosen  such  that  only  the  fundamental  transverse 
mode  oscillates,  and  thus  the  radiation  is  emitted  in  a  single 
transverse  lobe. 

If  we  limit  our  discussion  to  devices  that  have  koI  =  0  except 
for  the  active  layer  2,  then  (8)  assumes  a  form 

^ To  +  *0^00x^00]  g==0  (9) 

where  r2(y)  is  the  overlap  or  confinement  factor  of  the  mode 
power  to  the  active  layer.  A  similar  equation  was  originally 
introduced  by  Paoli  (21].  If  ne0  *  -iy0lk0  is  defined  as  the 
effective  index  in  the  absence  of  lateral  index  variations,  then 
the  net  lateral  effective  index  is 

n]  *  0'To/*o)*  +  rjCvKiO)-  (10) 

In  double-heterojunction  devices  ra(y)  may  be  as  low  as  0.3 
for  narrow  active  layers,  however,  for  LOC  lasers  T2  might  be 
less  than  0.1. 


III.  Effective  Index  Variations  of 
DHand  LOC  Structures 

In  this  section  the  lateral  variations  of  the  effective  index 
will  be  discussed  for  a  variety  of  nonplanar  DH  lasers.  These 
calculations  have  been  reported  earlier  for  DH  (22] ,  [23]  and 
LOC  [24]  devices  with  convex  curvature.  We  will  give  tome 
approximations  for  DH  lasers  and  extend  the  results  to 
LOC  devices  with  concave  curvature  of  the  heterojunction 
boundaries. 

For  simple  DH  structures  with  an  active-layer  thickness  d 
(Fig.  1),  the  effective  index  in  terms  of  the  normalized  width 


D-M(«j-«i),/*  is  [25] 

/!*#  »  bnj  +  (1  -  b)nf 

01) 

where 

,  .  D1  (y/974D*~  IV 

4  \  2  +  0*  / 

(12a) 

Fig.  3.  The  effective  index  variation  in  die  lateral  direction  for  both 
index  and  index  antiguides.  The  active  layer  varies  from  0.25  iitn 
to  0.1  i<m. 


has  been  approximated  from  the  confinement  factor  [26] 

r(D)-^  (12b) 

and  by  using  the  relationship  between  b,  T,  and  D  [27],  The 
approximation  formula  for  b  is  accurate  within  3  percent  over 
the  interval  0  <  D  <  2.5,  which  covers  all  cases  of  practical 
interest  in  AlGaAs/GaAs  and  InGaAsP/InP  devices.  Note  that 
for  InGaAsP/InP  structures,  where  D  is  virtually  independent 
of  X  [28]  (i.e.,  Dat6.6d),  one  obtains  an  approximation 
formula  for  b,  which  is  a  function  of  active-layer  thickness 
alone. 

In  Fig.  3  we  show  the  variation  of  the  effective  index  nt0  as 
a  function  of  y  for  both  convex  (index  guided)  and  concave 
(index-antiguided)  AlGaAs/GaAs  structures.  In  the  convex 
structure,  the  waveguide  width  varies  from  0.25  pm  to  0.1  Mm, 
while  the  concave  structure  is  of  opposite  variation;  the  refrac¬ 
tive  indexes  of  the  various  layers  are  n,  *3.4  and  n3  *3.6 
We  have  assumed  that  the  thickness  variation  follows  d  =  dt* 
(de  -  d0)  exp  (-y2lyl)  withy#  -  10pm. 

We  now  consider  the  waveguide  parameters  appropriate  to 
the  design  of  LOC  structures.  In  Fig.  4  the  confinement 
factors  corresponding  to  the  fundamental  and  the  first-order 
transverse  modes  are  shown  for  planar  LOC  structures  with 
d3  *  1  pm.  The  confinement  factors  of  these  transverse  modes 
(i.e.,  modes  in  the  plane  perpendicular  to  the  junction)  play 
roles  in  transverse  mode  selection  as  well  as  in  the  determina¬ 
tion  of  the  lateral  effective  index  variations  as  indicated  in  (9). 
At  y  *  0,  the  confinement  factor  of  the  fundamental  trans¬ 
verse  mode  must  be  larger  than  that  of  the  first-order  trans¬ 
verse  mode  if  the  device  is  to  operate  in  the  fundamental 
mode.  For  example,  if  d 3  *  0.2  pin  the  A1  concentration  of 
layer  3  must  be  larger  than  10  percent,  since  at  the  10  percent 
point  r~0.2  for  both  transverse  modes.  The  refractive 
indexes  of  the  AlxGat.xAs  layers  are  assumed  to  follow 
/r,  ■  3.6  -  0.62. r. 

A  peculiarity  of  LOC  structures  is  that  transverse  modes 
have  lateral  indexes  with  different  functional  dependencies. 
We  will  consider  specifically  a  LOC  device  which  is  grown  as  a 
concave-lens-like  structure.  Fig.  5  shows  the  various  regions 
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Fig.  4.  Confinement  factor  r  of  optical  power  to  the  active  region  for 
the  fundamental  transverse  mode  (solid  lines)  and  the  fust-order 
transverse  mode  (dashed  lines)  in  LOC  lasers  as  a  function  of  A1  con¬ 
centration  in  region  3  of  Fig.  2  (i.e.,  the  guide  layer).  The  graphs  are 
drawn  for  a  guide-layer  thickness  of  1  *im. 


LATERAL  DISPLACEMENT,  j  (MICRONS) 


Fig.  S.  A  cross  section  of  a  typical  nonplanar  (CDH-LOC)  laser.  The 
active-layer  thickness  decreases  with  lateral  position,  whereas  the 
guide-layer  thickness  increases.  The  effective  index  must  reflect  these 
opposite  dependences  on  lateral  displacement. 

with  thickness  variations  as  a  function  of  position,  as  discussed 
in  the  next  section. 

It  should  be  noted  that  the  effects  of  regions  2  and  3  on  the 
effective  index  oppose  each  other,  in  the  sense  that  one  pro¬ 
duces  focusing  while  the  other  defocusing.  As  mentioned 
previously,  the  A1  content  of  layer  3  affects  lateral  mode 
operation.  In  Fig.  6  we  show  the  lateral  index  variations  as  a 
function  of  y  for  the  fundamental  transverse  mode.  In 
Fig.  6(a),  where  x  *  9.7  percent,  the  lateral  wave  is  index- 
antiguided  when  gain  and  losses  are  not  considered.  On  the 
other  hand,  when  x  ■  11.3  percent,  the  wave  will  be  totally 
index  guided.  In  going  from  concavity  ( x  <  9.7  percent)  to 
convexity  (x>  11.3  percent)  the  intermediate  x  values  pro¬ 
duce  the  so-called  h'-shaped  waveguide  [29].  In  Fig.  7  we 
show  the  lateral  index  variations  as  seen  by  the  first-order 
transverse  mode;  note  that  this  mode  remains  index-antiguided 
for  all  A1  values  listed.  For  x  *  11.3  percent  the  fundamental 


transverse  mode  is  laterally  index-guided,  while  the  first-order 
transverse  mode  in  index-antiguided.  Similar  lateral  index 
variations  can  be  obtained  while  varying  the  layers’  thicknesses 
at  a  fixed  A1  concentration  [IS],  [29].  For  instance,  by  in¬ 
creasing  the  active-layer  thickness  to  0.25-0.3  um,  CDH-LOC 
devices  have  been  shown  [15],  [29]  to  become  totally  index- 
guided  for  the  fundamental  transverse  mode. 

IV.  Lateral  Wave  Confining  Structures  Grown 
on  Nonplanar  Substrates 

As  discussed  in  the  previous  section,  confinement  of  the 
optical  mode  can  be  realized  by  inducing  thickness  variations 
of  the  waveguiding  layers  along  the  plane  of  the  junction. 
Such  local  variations  in  thickness  can  be  generally  achieved 
.  while  depositing  material  by  liquid-phase  epitaxy  (LPE)  over 
nonplanar  substrates  [1]-[15].  We  show  in  Fig.  8  three  types 
of  laser  structures  that  can  be  obtained  by  LPE  over  topo¬ 
graphical  features  etched  into  the  substrate.  The  grown  layers 
have  variations  in  thickness  due  to  two  reasons:  1)  a  strong 
dependence  of  LPE  growth  on  local  surface  curvature  [2], 
[3],  [30],  [31],  and  2)  LPE-growth  sensitivity  to  the  degree  of 
substrate  misorientation  [3],  [30],  [31]. 

Fig.  8  displays  a  DH-type  device:  the  “ridge-guide”  CDH 
laser  [13],  [31]  [Fig.  8(a)],  and  two  LOC-type  devices:  the 
nonplanar  LOC  la«r  [10]  [Fig.  8(b)]  and  the  CDH-LOC  laser 
[11]  [Fig.  8(c)] .  In  the  DH  device,  as  in  most  DH  nonplanar- 
substrate  devices  [ lj— [4],  [6]-[8],  [12]-[14],  [32],  the 
active-layer  thickness  variation  solely  determines  lateral  wave 
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Fig.  7.  The  effective  index  as  a  function  of  lateral  position  for  the 
first  high-order  transverse  mode  using  the  structure  of  Fig.  S.  The 
waves  are  index-antiguided  for  all  A1  values  considered. 


confinement.  By  contrast,  in  LOC-type  devices  both  the  thick¬ 
ness  variations  of  the  active  and  guide  layer  play  a  role  in  mode 
confinement  and  selection.  The  nonplanar  LOC  structure 
[10]  has  convex-lens-shaped  active  and  guide  layers  as  a  result 
of  growth  over  a  channel.  A  more  complex  situation  is  the 
CDH-LOC  structure  for  which  growth  is  performed  above  the 
mesa  separating  a  pair  of  channels.  Growth  above  the  mesa  is 
slow  compared  to  growth  in  the  planar  areas  outside  the 
channels.  Thus,  after  the  growth  of  the  n-AlGaAs  confinement 
layer,  one  can  obtain  a  shallow  flat-bottomed  channel.  Then 
the  guide  layer  assumes  a  concave-lens-like  shape,  and  the 
subsequently  grown  active  layer  will  assume  a  convex-lens-like 
shape. 

It  should  be  stressed  that  uniform  current  flow  across  the 
lasing  area  is  “built-in”  for  CDH-type  devices  [Fig.  8(a) 
and  (c)]  since  in  such  structures  the  lasing  cavity  is  grown  on 
the  least  resistive  electrical  path  between  the  metallic  stripe 
contact  and  the  highly  conductive  substrate  [31].  By  contrast, 
in  structures  grown  over  channels  [l]-[5],  [9],  [10],  [32] 
(e.g.,  Fig.  8(b)] ,  the  current  has  the  tendency  to  be  focused 
towards  the  channel  shoulders,  and  thus  tight  current  confine¬ 
ment  to  the  lasing  area  (e.g.,  very-narTow-stripes,  preferential 
Zn-difTusion,  back-biased  junction)  must  be  provided  in  order 
for  the  fundamental  mode  to  be  excited. 


V.  Lateral  Mode  Behavior 

The  analysis  of  the  lateral  mode  structure  has  been  accom¬ 
plished  by  considering  several  types  of  dielectric  profiles  [33], 
[34].  These  ad  hoc  profiles  seldom  describe  the  actual  ones, 
but  an  understanding  of  the  real  modes  can  be  realized  by 
investigating  the  approximate  models,  in  each  of  the  lateral- 
index  functions  investigated  above,  it  should  be  noted  that 
with  the  assumed  thickness  variations  the  index  is  a  smooth 
function  of  lateral  position.  Of  all  analytical  profiles  which 
have  known  field  solutions,  the  hyperbolic  cosine  variation 
most  closely  matches  our  requirements  of  a  smooth  index 
profile.  The  index  profile  is  then  represented  by  the  following 
variation  of  the  dielectric  constant  x : 


kOO  =  k*  + 


k„-  Kb 


“*’(£) 


(13) 


where  xfl  -  nj  =  x(0)  and  nb=nb  =  x  (<*>);  na,  nb  are  refractive 
index  values  at  the  center  and  edges  of  the  distribution, 
respectively,  and  y0  is  a  lateral  displacement  characteristic  of 
the  index  lateral  variation.  (The  values  xa  and  Kb  may  be 
complex.)  It  is  interesting  to  note  when  xa  and  Kb  are  real 
that  if  Ax  =  Ka  -  Kb  >  0,  proper  modes  (trapped)  exist,  but  if 
Ax  <  0,  the  fundamental  mode  as  well  as  high-order  modes 
are  leaky.  Our  primary  interest  here  is  to  understand  both 
trapped  and  leaky  modes  in  concave  structures,  as  these  de¬ 
vices  have  excellent  power  output  and  spectral  character. 
Assuming  that  ku2  -  0  and  substituting  (13)  into  (9),  we  have 


dx* 


+  Ujjxj,  + 


fcjAx 


cosh2(y/y0) 


+  TJ  U  =  0 


(14) 


where  x  (>>)  is  the  effective  dielectric  constant.  The  solution  to 
(14)  is 


s  =  0, 1, 2  •  ■  •  (15) 

where  b,  =  b0  -  s  and  b0  =  (k\y% Ax  +  £),/J  -  ^  and  Cf(f) 
are  Gegenbauer  polynomials,  the  first  two  of  which  are 

<2( f)-  1  (16a) 

and 

Cf(f)«2Xf.  (16b) 


The  real  value  of  b,  ■  b0  -  s  determines  whether  the  s  mode  is 
proper  or  improper  (leaky).  Putting  b,*b',  +  ib',',  the  bound¬ 
ary  that  separates  proper  and  improper  modes  occurs  when 
b',  *  0.  From  (15)  it  is  seen  that  if  b',<0  the  mode  field 
strength  increases  without  limit  at  large  lateral  distances, 
whereas  b\> 0  gives  |g,(y)l-» 0  as  y  •* •».  The  boundary  de¬ 
fined  by  b't  *  0  is  satisfied  by  [35] 


A,/A«  ,  16s(s  +  l)(2s  +  1)*  -  (2kbybA"lt2)* 
4(2r  +  1)*  (2kbyb&!nl2)2 


(17) 


where  Ax  ■  A'  +  rA".  A  similar  result  has  been  derived  for 
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(c) 

F]g.  8.  (a)  Nonplanar  DH  User  structure  obtained  by  one-step  liquid -phase-epitaxy:  “ridge-guide"  constricted  double¬ 
heterojunction  (CDH)  later  [13],  [31],  nonplanar  LOC  later  structures  obtained  by  one-step  liquid-phase  epitaxy, 
(b)  nonplanar  large-optical  cavity  (NP-LOC)  laser  [10],  (c)  constricted-double-heterojunctkm  large-optical-cavity 
(CDH-LOC)  laser  [11]. 


step  index  variations  (36].  Fig.  9  shows  regions  where  the 
lateral  modes  are  proper  or  leaky  in  terms  of  the  complex 
dielectric  step  Ax  and  the  width  2y0.  The  propagation 
constant  is 

7*  *  -*o  [*6  +  hJ/^o-Po)’]-  (18) 

If  we  neglect  waveguide  losses,  index-antiguided  modes  have 
negative  Ax  values.  The  magnitude  of  the  negative  step  in¬ 
fluences  both  astigmatism  and  attenuation  of  the  mode.  Con¬ 
sider  the  fundamental  s*0  mode  when  -£  KklylAn  <0, 
b0  is  real,  and  thus  there  is  no  phase  variation  of  the  field 
along  the  facet.  Astigmatic  leaky  modes  occur  only  when 
klylA*  <  In  most  instances  the  attenuation  of  leaky 
modes  is  a  direct  reflection  on  power  loss  in  the  lateral  direc¬ 
tions  so  that  in  active  devices  a  mode  having  a  zero  net  at¬ 
tenuation  occurs  when  power  supplied  by  the  active  regions 
offsets  the  losses. 

Let  us  now  turn  to  waveguides  in  laser  structures  which  have 
gain/losses.  The  lateral  profile  must  be  modified  to  include 


the  dielectric  variations  occurring  in  the  active  layer.  From 
(9)  the  profile  is 

*0)  *  KoOO  +  T* OOx.-aO').  (19) 

Neglecting  temperature  and  strain  effects 

*„i(0)“  tg(0)  fit  lk0  (20a) 

and 

Kvj(*)* (20b) 

where  it  is  assumed  that  g(0)  is  the  gain  at  y*0  and 
<*]  «  -g(«»)  is  the  absorption  constant  of  layer  2  at  large 
lateral  distances.  In  the  unpumped  active  layer,  the  absorption 
coefficient  at  ~  100  -  500  cm'1.  If  we  make  the  ad  hoc 
assumption  that  rJ0')xl),0')  hu  a  lateral  dependence  similar 
to  KoOO,  then  x  O')  can  be  approximated  to  exhibit  the 
hyperbolic  cosine  variation.  The  dielectric  step  becomes 

Ax  -  -  nj  +  int  (f,(0)g(0)  +  Tjf-) a, )/k0  (21) 
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Fig.  9.  Cutoff  curves  for  the  lateral  modes  as  a  function  of  the  complex 

dielectric  step  and  waveguide  "width" >>0.  In  the  cross-hatched  region 

the  fundamental  mode  is  leaky  or  improper. 

with 

A'  =  n*  -  rt\  (22 a) 

and 

A"  =  n2  [r,  (0)  *(0)  +  r2  («)  Oj  ]  /k0 .  (22b) 

In  (22b)  the  total  value  of  the  imaginary  dielectric  step  is 
the  sum  of  two  positive  terms  because  the  gain  is  opposite  in 
sign  to  the  absorption  so  that  the  effective  gain  changes  from 
r2  (0)  g(0)  to  -  r2  (~)  a2 .  In  some  devices  the  current  distribu¬ 
tion  might  tend  to  be  uniform  in  the  region  containing  the 
dielectric  waveguide  so  that  the  gain  in  the  active  layer  varies 
according  to  the  thickness  of  the  active  layer  d2.  If  the  current 
has  inconsequential  lateral  variations  along  the  lateral  direc¬ 
tions,  the  gain  g(cm"‘)  is  written  as  [37]  (neglecting  lateral 
carrier  diffusion) 

*«45(/,ff/</2)-  190 

where  *  T)tJ  and  t)i  is  the  internal  efficiency,  7(kA/cml)  is 
the  injection  current  density,  and  d2  is  dimensioned  in 
microns.  Now  the  imaginary  part  of  the  dielectric  step 
becomes 

A"  -  n2  [r2(0)g(0)  -  r2(«)g(~)]/*0.  (22c) 


the  channel-substrate-planar  laser  [S]  have  only  trapped  lateral 
modes  which  lase;  the  structure  is  designed  so  that  the  real 
part  of  the  lateral  dielectric  step  is  positive.  In  the  CDH  laser 
discussed  here,  the  structure  is  fabricated  with  heterojunctions 
forming  a  concave  waveguide  [6]  that  has  an  active  region 
thickness  of  d2  =  0.16  pm  at  the  lasing  spot  (y  *  0)  and 
d2  =  0.28  pm  at  large  lateral  distances,  defined  in  our  formula¬ 
tion  as  y  *  Using  the  A1  concentrations  in  the  various 
layers,  we  have  na  ■  3.4866  and  nb  *  3.5213;  calculations 
give  IX0)« 0.445,  r(«)“0.71,  A' *  -0.243,  and  A"* 
(4.S4Jttf  +  24)  X  10*4,  where  the  gain  distribution  is  esti¬ 
mated  from  the  effective  current  density  which  is  almost  uni¬ 
form  over  the  waveguide  vicinity.  Assuming  a  value  y0  =  6 
pm,  the  condition  for  a  proper  mode  A"  *  (-  A')l,ilk0y0 
(crossing  from  the  curve  bounding  the  cross-hatched  region  of 
Fig.  9)  gives  /eff  ~  20  kA/cm*  which  corresponds  to  gain 
values  g( 0)  ~  5000  cm-1 .  Thus,  it  is  not  possible  for  this 
particular  leaky  device  to  reach  the  condition  of  a  fundamental 
proper  mode.  On  the  other  hand,  using  (18),  the  leaky  mode 
propagating  in  the  guide  will  reach  threshold  (assuming  end 
losses  ~30  cm*1 )  at  g( 0)  ~  600  cm'1  and  g(°°)  ~  260  cm*1 ; 
at  these  threshold  gain  values  b0  =  -0.33  +  f  22  which  would 
correspond  to  a  very  astigmatic  mode;  b'0,  the  real  part  of  b0 
means  the  intensity  increases  without  limit  at  large  y  values. 
However,  experimental  results  showed  that  such  devices  lased 
with  a  very  narrow  spot  (~2  pm)  at  the  facet.  There  are  two 
possible  explanations  for  the  discrepancy  between  our  calcula¬ 
tions  and  the  experimental  observations:  1)  in  the  waveguide 
region  at  y  *  0  the  actual  guide  was  relatively  flat  so  that  local 
ohmic  heating  in  the  region  ly  I  <  2  pm  produced  an  index 
increment,  thus  introducing  a  new  waveguide  parameter  y0 
(the  resulting  structure  would  in  fact  be  a  K'-guide),  and  2)  the 
possibility  of  “almost”  total  internal  reflection  of  the  mode  at 
the  waveguide  facet.  The  problem  with  the  first  explanation  is 
that  the  near  field  spread  should  be  much  larger  than  2  pm 
and,  in  fact,  the  laser  should  behave  more  like  a  stripe  geome¬ 
try  device.  The  possibility  of  total  internal  reflection  is  more 
likely  because  of  the  large  b'o  values  which  are  associated  with 
the  ray  directions  in  the  active  layer.  In  particular,  at  large  y 
values 


g(y)~e 


me~b*yfy»  e~lb*  yty> 


(23) 


Devices  exhibiting  steps  of  the  form  (22c)  are  commonly 
found  in  constricted  doubie-heterojunction  structures  unless 
tight  current  confinement  is  incorporated  in  the  laser  design. 

VI.  Device  Calculations 

Waveguide  parameters  will  now  be  calculated  for  index- 
antiguided  structures  (CDH  and  CDH-LOC  lasers)  to  illustrate 
the  utility  of  the  model. 

A.  "Leaky"  CDH  User  (  6J 

It  is  well  known  that  stripe  geometry  devices  have  modes 
supported  by  index  antiwaveguides,  and  that  the  modes  are 
essentially  trapped  because  of  the  lateral  gain  distribution. 
The  modes  are  not  leaky,  due  to  the  fact  that  the  optical  field 
in  the  active  layer  is  confined  more  or  leu  to  the  region  below 
the  stripe  contact.  Other  planar  active-region  devices  such  as 


The  two  plane  wave  components  in  region  2  which  form  the 
optical  mode  with  respect  to  the  x  direction,  are  propagating 
at  an  angle  with  respect  to  the  z  axis.  Thus,  the  net  propaga¬ 
tion  vector  xfix  +yf}y  +?&  of  a  single  “plane  wave”  in  the 
active  layer  has  fix  «  *0(*i  _  «3rf(y)),/* ,  Py-b'o/yo.  «nd 
Pz  =  Im(>),  determined  from  (18).  Combining  the  above 
components,  rays  in  the  active  layer  region  with  d2  =  0.28  pm 
propagate  at  an  angle  (with  respect  to  the  z  axis)  0  ~  14.2°, 
which  is  near  the  critical  angle  for  total  internal  reflection  at 
the  facet.  Further,  if  these  rays  escape,  it  is  doubtful  that  they 
entered  the  microscope.  This  explains  why  no  radiation  leak¬ 
age  was  observed  for  the  ‘leaky -guide”  CDH  device  [6] .  Only 
radiation  propagating  along  the  axis  could  be  detected,  and 
that  provided  a  relatively  narrow,  centered  beam  (0(  *  20s). 
As  the  width  d2  drops  from  0.28  pm,  the  value  of  0X  increases 
and  the  corresponding  ray  angles  increase.  But  near  y  *  0  the 
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plane  waves  change  directions,  thus  increasing  the  possibility 
for  light  escaping  from  the  facet. 

B.  CDH-LOC  Lasers 

The  geometry  of  CDH- LOC  devices  is  given  in  Fig.  5;  we 
will  use  the  structure  with  9.7  percent  A1  in  layer  3.  This  is 
an  index-antiguided  structure  that  can  operate  in  a  high -order 
transverse  mode,  but  the  fundamental  is  assumed.  Aty  =  0, 
d2  =  0.2  um  and  d3  =  1.3  ftm,  whereas  at  large  y,  dt  =0.15 
Aim  and  d3  =  2.5  Atm.  The  confinement  factors  are  r2(0)  = 
0.157  and  r2(®®)  =  0.023.  The  refractive  indexes  na  =  3.5363 
and  nb  =  3.5373  are  used  to  get  A'  =  -0.007.  The  value 
y0  =  5.3  Atm  is  estimated  from  lateral  points  where  the  refrac¬ 
tive  index  changes  go  through  one  half  of  their  total  change, 
i.e.,  if  yH  represents  the  point  of  half  variation,  y0  -  0.759  y*. 
Assuming  a  uniform  current  distribution  over  the  lasing  por¬ 
tion  of  the  waveguide,  we  find  for  a  trapped  or  proper  mode 
to  exist,  J,tt  ~~  2.5  kA/cm3  which  corresponds  to  g(0)  ~  370 
cm'*  and  g(«®)  ~  560  cm"* .  Even  though  g(®®)  >g(0),  the 
effective  gains  are  modified  by  the  confinement  factors  of  the 
various  lateral  points.  At  Jt(f  =  2.5  kA/cm3 ,  the  fundamental 
mode  attenuation  coefficient  a0  ~  13  cm"1 ,  but  at  threshold, 
a0  =  -30  cm"1  and  Jtn  ~  3  kA/cm3 ,  which  give  g( 0)  =  485 
cm'1 ,  g(®®)  =  710  cm'1 ,  and  b0  =  0.16  +  /3.4.  At  these  small 
b'o  values,  the  rays  in  the  active  layer  are  determined  pre¬ 
dominantly  by  Px  and  02 . 

As  a  final  example,  we  consider  a  CDH-LOC  laser  [15] 
that  supports  index-guided  modes.  This  structure  has  n,  = 
3.418,  n2  =  3.6,  n3  =  3.502,  and  n4  =  3.44  and  d 2  ~  0.27  A«m, 
dj  ~  1.8  pm  at  y  -  0  and  d2  ~  0.2  Aim,  d 3  2.1  Aim  at  y  ~  10 

Aim.  Computations  give  na  =  3.5297,  nb  =  3.5153,  F(0)  = 
0.6157  and  I\°°)  =  0.4282.  From  geometrical  estimates, 
y*  ~  7  Aim  giving  y0  ~  5.3  Aim.  Again  assuming  a  uniform 
current  distribution,  threshold  occurs  when  g( 0)  ~  60  cm"1 
and  £(®°)  ~  150  cm'1.  At  these  gain  values,  b0  =  12.3  -  i0.076. 
Actually,  the  gain  distribution  in  this  structure  plays  little  or 
no  role  in  the  lateral  mode  shape.  Nevertheless,  the  large 
charge  density  differential  between  the  points  y  =  0  and  y  =  10 
Aim  will  cause  a  relatively  large  diffusion  current  directed  to¬ 
ward  y  =  0.  It  is  interesting  to  note  that  our  model  holds  for 
cases  such  as  when  the  gain  is  suppressed  at  y  =  0  and  increases 
with  y;  this  is  opposite  to  the  CDH  structure  discussed  above. 

VII.  Radiation  Patterns 

The  optical  radiation  pattern  can  be  determined  once  the 
near  field  distribution  is  known.  The  most  elementary  ap-. 
proach  and  the  one  given  here  is  to  calculate  the  Fourier 
transform  of  the  near  field  mode.  This  method  is  reasonably 
accurate  for  proper  modes  where  there  is  little  internal  mode 
conversion  at  the  facet.  In  the  case  of  the  “leaky”  CDH  de¬ 
vice  [6]  discussed  earlier,  portions  of  the  fundamental  mode 
in  the  regions  y  >y0  might  see  a  very  high  reflection  coeffi¬ 
cient  at  the  facet.  Thus,  little  light  would  escape  from  the 
device,  with  most  coming  from  the  region  around  y  =  0. 
Under  these  conditions,  much  of  the  reflected  light  would  be 
distributed  in  the  waveguide  radiation  modes  and  be  dissipated 
internally. 
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Fig.  10.  The  half-power  be  am  width  9.  as  a  function  of  the  waveguide 
width  2y„  for  the  CDH-LOC  laser  of  [  1 5 ) . 

Since  the  near  field  is  separated  into  /(*,  y)  and  g(y)  via  the 
effective  index  method,  the  far  field  can  be  similarly  divided 
into  independent  perpendicular  angle  dependencies.  We  con¬ 
centrate  here  on  the  lateral  pattern  F(0)  given  by 

F(8)  =  J  g(y)e‘xy  dy  (24) 

where  x  =  *o  sin  0  and  0  defmes  the  lateral  angle  from  the 
facet  normal.  (The  obliquity  factor  has  been  dropped.)  For 
the  sake  of  briefness,  we  consider  only  the  fundamental  mode 
pattern  in  the  lateral  direction.  Using  g(y)  =  cosh'*®  (y/y0) 
-»,-z 

HO)  *  Wo/2  +  «yoX/2)  IX*o/2  -  i>oX/2)  (25) 

where  IXw)  is  the  gamma  function  [not  the  confinement  factor 
m  (6),  (9),  (12a),  and  (19)].  Contour  and  relief  maps  [38]  of 
the  gamma  function  with  complex  arguments  can  be  used  to 
analyze  the  behavior  of  the  far  field  patterns.  The  gamma 
function  F(w)  has  a  pole  at  w  *  0.  If  we  put  w  =  u  +  iu  with 
u  -  constant  and  vary  v,  the  gamma  function  peaks  at  o  =  0 
and  varies  monotonically  from  0  =  0.  In  the  far  field  pattern 
of  the  fundamental  mode,  the  two  gamma  functions  actually 
form  different  peaks  at  symmetrical  positions  about  k0  sin 
0  =  0.  If  the  value  b'o  is  sufficiently  small  the  two  peaks  com¬ 
bine  to  form  a  single  peak  at  0  =  0.  But  in  the  case  when  the 
two  peaks  can  be  separated  and  b'o  is  large,  peaks  occur  at 

b" 

sin  0  =  ±  — — .  (26) 

k0y0 

As  the  value  of  b'o  -*•  0  from  the  positive  side,  the  two  peaks 
increase  without  limit.  The  value  of  b'o  is  dependent  upon  the 
dielectric  step  Ax  and  k0y0.  If  Re  [Ax]  is  positive  and  large, 
b'd  ~  0,  but  if  Re  [Ax]  is  large  and  negative,  then  for 
*oyo  »  0,  b'o  ~ koyofal  ~  ns)l/2  *o  that  far  field  peaks 
occur  at 

sin  0  =  ±(nl  -  n J)'/1  (27) 

which  allows  us  to  estimate  the  index  step  for  index-antiguided 
modes.  For  large  bo  values,  Stirling’s  approximation  for  the 
gamma  function  can  be  employed. 
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Fig.  11.  Theoretical  experimental  pattern  of  the  CDH-LOC  User  of 
[15].  The  optimum  fit  occurs  atyo  *  6.2S  pm. 

To  illustrate  the  use  of  (25),  we  calculate  the  radiation  pat- 
tern  of  the  index-guided  CDH-LOC  (15]  discussed  earlier. 
Using  the  index  calculation,  Fig.  10  shows  the  half-power 
beamwidth  of  the  fundamental  lateral  mode  as  a  function  of 
2 yn-  At  2 yn  ~  14  fan  we  find  that  0^  ~  8°,  which  is  in  good 
agreement  with  the  experimental  patterns.  In  Fig.  1 1  we  plot 
both  theory  and  experiment. 

VIII.  Conclusion 

The  effective  index  method  has  been  applied  to  both  CDH 
and  CDH-LOC  lasers  to  analyze  their  lateral  mode  behavior. 
The  LOC  devices  have  shown  that  lateral  mode  behavior  is 
strongly  influenced  by  transverse  mode  operation  when  the 
optical  cavity  contains  AlGaAs  layers  with  relatively  large 
aluminum  concentrations.  This  occurs  in  geometries  where  the 
optical  cavity  increases  in  thickness  laterally.  In  devices  with 
relatively  small  active  layer  thicknesses  all  transverse  modes 
are  index-antiguided  in  the  lateral  direction. 

Lateral  modes  are  analyzed  using  the  smoothly  varying 
hyperbolic  cosine  distribution.  The  effects  of  the  active 
region  gain  are  included  in  the  model  to  indicate  when  leaky 
waveguide  modes  become  proper. 
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Transverse-mode  selection  is  characterized  for  GaAs/AlGaAs  double  heterojunction  lasers  from 
optical  field  and  electron/hole  interaction.  The  electron/hole  distribution  determined  from  a 
solution  of  the  ambipolar  diffusion  equation  provides  the  necessary  information  about  gain/mode 
coupling  to  predict  the  current  at  threshold.  Lasing  power  out  versus  current  solutions  provide 
information  about  internal  differential  quantum  efficiency.  Theory  is  matched  to  experiment  for  a 
multimode  laser  with  one  heterojunction  having  a  very  small  index  step.  It  is  found  that  the  laser’s 
characteristics  over  a  temperature  and  current  range  are  predicted  by  adjusting  the  active-layer 
refractive  index  as  determined  from  far-field  measurements. 


PACS  numbers:  42.55.Px,  42.60.By 


I.  INTRODUCTION 

In  recent  years,  there  has  been  considerable  effort  spent 
fabricating  various  semiconductor  laser-device  geometries 
for  control  of  mode  operation.  In  GaAs/AlGaAs  devices, 
the  transverse-mode  operation  is  usually  governed  by  the 
growth  of  various  layer  thicknesses  and  index  steps  at  the 
grown  heterojunctions.  Double  heterojunction  lasers  which 
confine  the  optical  field  to  the  thin  active  layers  usually  oper¬ 
ate  in  the  fundamental  mode.  Furthermore,  the  lasing  mode 
is  independent  of  drive  because  small  index  changes  due  to 
injection  have  negligible  effects  on  mode  selection.  While 
transverse-mode  lasing  can  be  well  controlled  by  epitaxial 
heterojunctions,  lateral-mode  operation  in  contemporary  la¬ 
sers  is  drive  sensitive.  Stripe-contact  devices'-3  have  lateral 
guides  defined  only  by  the  current  density  which  affects  the 
gain  distribution  while  more  sophisticated  structures4-7 
have  grown  layers  whose  thicknesses  are  functions  of  lateral 
positions.  Thus,  in  these  devices  the  lateral  modes  are  shaped 
by  the  lateral  effective  index  as  well  as  gain  variations. 

Because  the  gain  distribution  plays  a  major  role  in  de¬ 
fining  mode  stability  in  lasers,  it  is  important  to  understand 
the  mechanisms  of  electron/hole  transport  in  active  layers. 
In  this  paper,  we  discuss  the  transverse-mode  operation  of 
lasers  which  are  affected  by  the  electron/hole  distribution. 
For  the  first  time  the  ambipolar  diffusion  equation  is  solved 
for  the  active  layer  in  an  injection  laser.  The  direct  conse¬ 
quence  of  this  solution  is  the  ability  to  predict  crossover  of 
competing  modes  in  the  active  layer.  In  contrast  to  ambipo¬ 
lar  diffusion,  assuming  simple  electron  injection  does  not 
explain  this  phenomena  in  sufficient  detail. 

In  this  work  we  investigate  mode  selection  from  the 
standpoint  of  the  exact  selection  process.  This  includes  mode 
gain  coupling  and  differential  quantum  efficiency.  We  con¬ 
centrate  on  the  transverse  (perpendicular  to  the  junction 
plane)  modes  in  a  double  heterojunction  laser.  The  excess 
carrier  distribution  is  assumed  ambipolar  and  a  complete 
solution  is  obtained.  From  mode/gain  coupling,  it  is  possible 


to  predict  the  threshold,  and  from  the  slope  efficiency,  it  is 
possible  to  predict  mode  crossover. 

The  laser  we  apply  this  to  is  a  broad-area  double  hetero¬ 
junction  GaAs  device  with  a  small  index  step  at  one  of  the 
heterojunctions.  The  laser  is  heated  over  a  temperature 
range  and  the  first  two  modes  are  investigated.  It  is  found 
that  cavity-mode  preference  is  very  sensitive  to  temperature 
and  current  level.  Application  of  the  diffusion  equation  is 
used  to  describe  the  crossover  of  the  first  two  modes  and  it  is 
established  that  the  mode  with  the  highest  slope  efficiency 
will  dominate.  Far-field  patterns  are  described  by  standard 
multilayer  waveguide  techniques.  Finally,  a  calculation  is 
made  to  estimate  the  effect  of  free  carriers  on  the  refractive 
index  in  the  material. 

II.  THEORY 

A.  Approach 

Excess  carrier  movement  in  the  active  layer  of  an  injec¬ 
tion  laser  is  described  using  the  ambipolar  diffusion  equa¬ 
tion.  Quasi-Fermi  level  continuity  as  well  as  electron  and 
hole  current  continuity  is  imposed  at  each  interface.  The 
quasi-Fermi  level  location  is  calculated  using  a  nonparabolic 
I*  i  c  extremum,  and  parabolic  XIC  and  f': uv  extrema.  The 
boundary  condition  is  set  up  in  such  a  way  as  to  eliminate  the 
need  for  a  solution  within  the  space  charge  layer  surround¬ 
ing  each  interface.  The  spontaneous  recombination  lifetime 
is  calculated  without  imposing  quasimomentum  conserva¬ 
tion  for  parabolic  bands.  Above  threshold,  a  stimulated  re¬ 
combination  term  is  included  in  the  ambipolar  diffusion 
equation  and  the  gain  is  assumed  to  be  pinned  at  the  thresh¬ 
old  value.  The  lasing  modal  shape  for  the  transverse  modes 
(perpendicular  to  the  metallurgical  junctions)  is  calculated 
by  solving  Maxwell’s  equations  for  a  multilayer  dielectric 
wave  guide. 

B.  The  excess  carrier  distribution 

The  diffusion  equation  evolves  along  the  standard 
route.  We  demand,  as  is  customary,  that  each  carrier  type 
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satisfies  the  continuity  equation8 

Dp  V2P  -  ?•( pE  i  +  gp  -  p/rp  =  dp/dt,  ( la) 

£>„ V2n  +  p„  V-(nE )  +  g„  -  n/r„  =  dn/dt,  (lb) 

where  n,  p  =  total  electron,  hole  carrier  density;  D„,  Dp 
=  electron,  hole  diffusion  constant;  g„,gp  =  electron,  hole 
generation  rate;  and  r„ ,  =  electron,  hole  mean  recombin¬ 

ation  time.  We  further  demand  that  Poisson's  equation  be 
satisfied: 

V-E  =  (4  irp/e),  (2) 

where  p  is  the  space  charge  density  and  e  is  the  dielectric 
constant. 

Rather  than  solve  Eqs.  (1)  and  (2)  exactly,  we  take  an 
approximate  route  where  we  assume  that  the  internal  field 
between  the  electron  and  hole  charge  densities  is  strong 
enough  to  guarantee  charge  neutrality.  Although  this  condi¬ 
tion  is  violated  in  the  space  charge  layer  surrounding  the 
heterojunctions,  we  can  formulate  the  problem  in  such  a  way 
that  only  diffusion  processes  outside  the  heterojunctions 
need  be  considered.  Under  this  assumption,  the  excess  elec¬ 
tron  and  hole  densities,  Sn  and  Sp,  respectively,  must  be 
equal: 

Sn  =  n~n0=p—p0-Sp.  (3) 

Here,  /»„,  p0  designate  the  thermal-equilibrium  carrier  densi¬ 
ties.  If  we  take  \dn0)/[dx)  —  {dp0)/{dx)  =  0,  we  can  substitute 
Eq.  (3)  into  Eqjl)  to  obtain  two  equations  for  the  two  unk¬ 
nowns  6p  and  E.  The  terms  involving  V-E  can  be  eliminated 
by  multiplying  Eq.  (la)  by  np,  and  Eq.  (lb)  by ppp  and  add¬ 
ing  to  obtain9 

D  *V2{6p)  -  p*E-V{6p)  4-  g'  —  —  =  ,  (4) 

t  at 

where 


Sp  =  excess  carrier  density, 

D  .  _  +PHpDp 

nPn  +  PPp 

HUn  +  PPp 

Sp  _Po  +  &p  Pq  n0  +  Sn  n0 


g'„  &'p  =  excess  carrier  generation  rates,  and 
r  =  excess  carrier  recombination  lifetime. 


We  have  implicitly  assumed  that  the  net  recombination  rate 
(8p)/r  is  the  same  for  electrons  as  for  holes,  as  it  must  be.  We 
have  further  assumed  that  the  equilibrium  generation  rate  is 
equal  to  the  equilibrium  recombination  rate  for  both  elec¬ 
trons  and  holes. 

For  the  ambipolar  diffusion  Eq.  (4),  observe  that  the 
ambipolar  diffusion  rate  for  excess  carriers  is  determined  by 
the  average  value  of  D„  and  Dp ,  weighted  by  np„  and  ppp .  In 
other  words,  the  two  oppositely  charged  clouds  interact  with 
each  other  through  the  internal  field  in  such  a  way  as  to 
satisfy  Poisson’s  equation  and  must  then  diffuse  at  a  natural¬ 
ly  acceptable  rate.  This  is  crucial  because  of  the  very  large 
injected  charge  densities  encountered  in  injection  lasers 


(~2x  1018  cm-3). 

The  electric  field  in  Eq.  (4)  refers  to  both  the  internal 
and  applied  fields.  Under  the  conditions  encountered  in  an 
injection  laser,  these  fields  exert  only  a  small  direct  influence 
on  carrier  distribution,  and  as  such  can  be  ignored  in  Eq.  (4). 
In  other  words,  the  electric  field  does  not  alter  the  actual 
excess  carrier  spatial  distribution,  but  rather  only  causes  car¬ 
riers  to  drift.  This  is  the  case  if  the  drift  length  is  much  less 
than  the  diffusion  length,  i.e., 

P*Et<(D*t)u1.  (5) 

This  condition  is  met  in  the  solutions  ultimately  obtained. 
The  major  effect  of  the  internal  field,  that  of  charge  neutra¬ 
lity,  is  still  accounted. 

If  we  assume  the  only  excess  carrier  drain,  other  than 
spontaneous  recombination  Sp/r,  is  stimulated  recombina¬ 
tion,  then  for  an  infinitely  long  cavity  and  time-averaged 
values10 


r_  m\  i  gp(x) 

dt  i f,  hv 


(6) 


where  P  (x)  =  time-averaged  power  density,  g  =  gain,  rj, 
—  stimulated  quantum  efficiency,  and  hv  =  photon  energy. 
For  100%  carrier  confinement,  rj,  is  the  internal  differential 
quantum  efficiency.  In  these  calculations  we  put  if,  =  1. 
Equation  (4)  becomes 


D*V2(Sp)-(Sp)/r  =  r, 


(7) 


where  2)  •  =  ambipolar  diffusion  constant,  r  =  spontaneous 
recombination  lifetime,  and  r  —  stimulated  recombination 
rate.  The  value  r  is  equated  to  zero  in  the  region  where  g  <  0. 
The  spontaneous  recombination  lifetime  is 


r  =  l/[B(n0  +  p0  +  Sp)],  (8) 

where  B  is  1.3x  lO^'W/sec."-'2 


C.  The  boundary  conditions 

At  each  boundary  within  the  device  there  exists  a  space 
charge  layer  with  a  width  severely  reduced  under  conditions 
of  forward  bias.  Because  quasicharge  neutrality  is  violated 
within  this  layer,  we  match  across  the  layer  under  the  as¬ 
sumptions  it  is  is  thin  enough  and  internal  fields  are  strong 
enough  that  carriers  are  swept  through  without  significant 
recombination.  Hence,  we  demand  that  the  electron  and 
hole  current  densities  be  the  same  on  each  side  of  the  space 
charge  layer.  These  densities  are  written 

7.  =  qnpHE  +  qD„  V(6p),  (9a) 

~jp  =  <U>HpE  ~  qDf  V{6p).  (9b) 

Here,  £_is  the  total  field  and  q  is  the  electron  charge.  Al¬ 
though  E  can  be  ignored  in  Eq.  (7),  where  it  has  a  small  effect, 
it  must  be  included  in  the  boundary  conditions  because  both 
drift  and  diffusion  terms  are  significant. 

In  addition  to  continuity  of  current  densities,  we  also 
impose  continuity  of  hole  and  electron  quasi-Fermi  levels 
across  each  layer  interface.  Although  it  is  clear  this  align¬ 
ment  of  Fermi  levels  implies  band  bending,  i.e.,  space 
charge,  we  assume  an  abrupt  heterojunction  and  ignore  such 
effects. 
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0.  The  electromagnetic  field  distribution 

The  standard  multilayer  slab  dielectric  waveguide 
model  is  used  to  describe  the  transverse  far-field  pat¬ 
terns.  ,3  M  Standard  calculations  yield  propagation  constants 
y  —  a  -+■  j(3  and  fields.  For  instance,  for  a  TE  mode,  we  can 
express  the  mth-mode  field  distribution  for  the  /th  layer  as 

E‘Jx)  =  Em  cos(h'„x  +  d'J,  (10) 

where  h  ‘m  is  the  complex  eigenvalue  and  x  is  the  coordinate 
across  the  active  region.  The  overall  field  solution  is  obtained 
by  matching  the  field  components  at  each  interface.  The  far- 
field  pattern  then  becomes  proportional  to  the  Fourier  trans¬ 
form  of  the  near  field  on  a  lasing  facet. 


E.  The  diffusion-equation  solution 

Assume  a  three-layer  device  with  the  central  layer  ac¬ 
tive  (see  Fig.  3).  Also,  assume  the  mode  is  propagating  down 
the  z  axis  with  x  perpendicular  to  the  guide.  The  excess  car¬ 
rier  diffusion  equations  for  this  case  become 

0,  ;=  1.3.  (lla) 

dx1  \L  f)2 


i35eL— &— i=  2,  (lib) 

dx 2  {L?)2  vM 

where  L  *  =  (D  *Y)l/2,g(x)  is  the  localized  gain  and  Pm  |x)  is 
the  optical  Poynting  vector  component  along  the  propaga¬ 
tion  direction. 

These  equations  are  linearized  by  replacing  gfx)  with  its 
weighted  average  within  the  active  region10 

S1mg\x)E2Jx)dx  gmNm 
g  SE2m(x)dx  Tm 

where  Nm  is  the  field  normalization  constant  and  Tm  is  the 
mode  intensity  integrated  over  the  active  layer.  It  should  be 
noted  that  g(x)  =  0  outside  the  active  region.  Further,  as¬ 
suming  no  losses  in  layers  external  to  the  active  region,  g„ 
represents  the  propagation  gain  of  the  mth  mode.  With  D  * 
and  r  assumed  to  be  piecewise  constant  layer-to-layer,  Eq. 
(lib)  becomes 


d2(Sp) 

dx2 


cos 2(h  lx  +  <t>: 


where 


PJx)dx. 


Ignoring  the  minor  effect  of  the  layer  3-4  interface  on  the 
carrier  distribution,  the  solution  to  the  diffusion  equation  in 
the  various  layers  becomes 

Sp3{x)  =  G3  cosh(x/I  ?)  4-  F3  sinh(x/L  ?),  (14a) 

Sp2  =  Gj  cosh(x/L  J)  -I-  F2  sinh(x/L  J) 


6pt(x)  =  G,  cosh(x/L  f )  +  F,  sinh(x/£  f ).  . 


With  these  eigenfunctions,  the  solution  follows  directly 
from  matching  the  current  density  and  quasi-Fermi  levels  at 
each  interface.  The  Fermi-level  continuity  condition  is  ap¬ 
plied  using  Kane’s  nonparabolicity  for  the  riC  minimum 
and  a  parabolic  expression  for  both  riiV  and  the XIC  mini¬ 
ma.  The  effect  of  aluminum  concentration  on  carrier  mobil¬ 
ity,  effective  mass,  and  band  gap  is  included.  The  volume 
recombination  rate  is  related  to  the  optical  gain  of  the  laser  in 
linear  fashion  using  coefficients  obtained  by  integrating  over 
parabolic  bands  without  quasimomentum  conservation.  The 
recombination  lifetime  is  found  in  a  similar  manner.  The 
diffusion  coefficients  are  calculated  using  the  generalized 
Einstein  equation.  The  complete  formulation  of  the  problem 
is  presented  elsewhere.13  The  necessary  roots  are  calculated 
on  a  computer  using  a  real-root  searching  routine.  The  final 
solution  provides  current-voltage  curves,  power  output-cur¬ 
rent  curves,  threshold  current  density  of  each  mode,  and 
finally  electron  and  hole  injection  efficiencies  into  the  active 
layer. 

F.  Validation  of  numerical  results 

This  theory  has  been  applied  to  double  heterojunction 
lasers,  single  heterojunction  lasers,  and  separate  electron- 
optical  confinement  lasers.  Close  agreement  has  been  found 
with  experimental  values.  Two  typical  results  for  a  double 
heterojunction  laser  are  shown  in  Figs.  1  and  2.  The  first 
shows  threshold  current  density  as  the  active  region  thick¬ 
ness  is  reduced,  and  the  second  shows  threshold  as  a  function 
of  the  aluminum  concentration  in  the  passive  regions.  The 
slope  of  Fig.  1  is  S.2  kA/cm2/xm,  a  value  which  compares 
closely  with  available  data.1617 


CAVITY  WIDTH  (MICRONS) 

FIG.  1.  Calculated  bole  injection  efficiency,  relative  peak  power,  and 
threshold  current  density  of  a  symmetric  DH  laser  as  a  function  of  cavity 
width.  The  active  layer  doping  isp- 10”  cm  " 1  and  the  passive  n-layer  doping 
is  10"cm-,.Thepassive^layerhasaSxl0"cm~>dopinglevd.TbeAlAs 
percentage  in  the  outer  layers  is  20%. 
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AIAs  (PERCENT) 

FIG.  2.  Calculated  threshold  current  density  as  a  function  of  AlAs  percen¬ 
tage  in  the  GaAlAs  passive  p  layer  of  a  symmetric  DH  laser  with  a  l^im 
wide  GaAs  (p-10n  cm-1)  active  layer.  Three  passive  p-layer  doping  levels 
are  shown.  The  pasrive  n-layer  doping  level  is  1 X 10"  cm-1. 


III.  EXPERIMENT 

The  lasers  characterized  in  this  experiment  are  broad- 
area  double  heterojunction  lasers  grown  by  liquid  phase  epi¬ 
taxy  methods.  Figure  3  gives  the  geometry  of  one  such  laser, 
PL174-1-31.  Other  lasers  of  this  type  were  characterized, 
but  PL174- 1-31  was  selected  for  its  transverse  modal  behav¬ 
ior  with  temperature  due  to  the  small  refractive  index  step  at 
one  heterojunction.  Estimating  the  refractive  index  of  the 
aluminum  layers  is  accomplished  using  experimental 
data1*-19  at  a  lasing  wavelength  A  =  9000  k.  This  gives 
4rj  =  0.62x,  where  x  is  the  fraction  of  aluminum  in  the  solid. 
The  laser  has  a  cross  section  of  214x356  fim,  where  the 
latter  is  the  cavity  length.  The  chip  thickness  is  77  fim.  The 


TRANSVERSE  DIRECTION  x 


FIG.  3.  Geometry  of  PL174-1-3I.  For  the  excess  carrier  calculation,  only 
layers  1, 2,  and  3  are  used. 
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laser  is  mounted  p-side  up  on  the  header  and  has  no  reflective 
coating  on  the  facets. 

The  device  is  driven  with  250-nsec  pulses  with  a  0.01% 
duty  cycle  to  minimize  internal  heating.  While  driven,  the 
laser  is  heated  over  temperature  range  of 20-40  *C.  The  tem¬ 
perature/current-drive  characteristics  of  the  laser  are  char¬ 
acterized  by  its  far-  and  near-field  dependence. 

To  measure  the  far  field,  the  laser  is  placed  on  a  rotating 
platter  with  an  axis  turned  by  a  2-rph  timing  motor.  The 
laser’s  output  is  collected  by  monochromator  and  amplified 
by  an  RCA  7102  photomultiplier  tube  (PMT).  Typical  mea¬ 
surement  parameters  are  2-A  slit  width  and  750-V  bias  for 
the  PMT.  The  output  ofthe  PMT  is  fed  to  a  lock-in  amplifier 
triggered  by  a  reference  signal  from  the  pulser.  An  HP  7560 
A  log  converter  and  HP  2470  AX  10  dc  multiplier  process 
the  signal  for  the  y  axis  of  a  chart  recorder.  An  impedance 
matching  network  allows  they  axis  to  be  calibrated  to  three 
decibels  per  inch  deflection.  The  calibration  is  good  for  over 
five  decades  of  response.  The  speed  of  the  x  axis,  combined 
with  the  2-rph  timing  motor  presents  24*  rotation  per  inch 
deflection  for  the  x  axis.  An  alternative  presentation  is  to 
load  the  PMT  with  a  50-/2  terminator  and  display  the  output 
pulse  on  a  sampling  oscilloscope  with  a  50-12  input.  This 
provides  a  time  decomposition  of  the  light-output  pulse  for 
ready  comparison  to  the  current-input  pulse.  Patterns  were 
typically  made  with  J  —  l.l/th . 

The  measurement  scheme  for  the  near  field  is  accom¬ 
plished  by  attaching  the  laser  mount  to  a  micropositioner 
which  has  a  three-dimensional  adjustment.  A  microscope 
objective  with  160  magnification  is  mounted  opposite  the 
laser  facet  and  a  video  camera  collects  the  output  ofthe  ob¬ 
jective  lens  and  displays  the  image  on  a  monitor. 

Peak  threshold  current  values  for  laser  PL174-1-31  are 
shown  in  Fig.  4.  Triangles  denote  the  fundamental  mode 
whereas  crosses  mark  the  second.  In  the  heat-sink  tempera¬ 
ture  range  20-28  *C,  only  the  second  mode  is  present;  for  28- 
42  *C,  both  fundamental  and  second  modes  are  present;  and 
above  42  *C,  only  the  fundamental  mode  propagates.  The 
lasing  wavelengths  for  the  two  modes  are  shown  in  Fig.  5. 

The  multilayer  waveguide  model  is  used  to  describe  the 
transverse  far-field  patterns  of  the  modes.  Figure  6  is  an  ex¬ 
ample  of  a  match  of  the  second  mode  at  20  *C.  The  solid 
curve  is  the  measured  pattern,  and  the  dashed  curve  is  the 
match.  The  bump  on  one  shoulder  at  21.6*  of  Fig.  6  is  the 
result  of  light  coupling  to  the  lossy  substrate.  This  bump  can 
be  accurately  positioned  to  match  the  experimental  pattern 
by  applying  Snell’s  law  at  the  substrate-air  interlace.  The 
internal  angle  from  the  normal  to  the  facet  O'  is  given  by 
tan  0  '~0  '~/i  *in  0.  Here,  h  « is  the  imaginary  part 

of  h4,  the  real  propagation  constant  /?  =  Im{ y\ ,  0  is  the  ex¬ 
ternal  angle  to  the  normal  to  the  facet,  and  ij4  is  the  refractive 
index  of  the  substrate.  The  relative  intensity  or  amplitude  of 
the  bump  is  adjusted  by  the  separating  width  dy 

0"  is  defined  as  the  angle  separating  the  two  major 
lobes  of  the  second  mode.  The  solid  line  in  Fig.  7  shows  how 
6"  changes  over  the  temperature  range.  The  dashed  line  in 
Fig.  7  is  the  multilayer  0„  versus  the  active-region  refractive 
index  rft.  Over  the  entire  range  20-42  *C  the  active-region 
index  depression  necessary  to  match  experimental  evidence 
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FIG.  4.  Peak  threshold  current  vs  heat-sink  temperature.  The  triangles  are 
data  points  for  the  fundamental  mode  and  the  crosses  for  the  second  mode 
of  PL174-1-31. 


ANGLE  FROM  NORMAL  (DEG  ) 

FIG.  6.  Experimental  (solid  line)  and  theoretical  (dashed  line)  far-field  pat¬ 
terns  for  the  seoond  mode  of  PL174-1-3I  at  20 f. 


is  6ij  =  —  0.002.  The  fundamental  mode  also  is  marked  by 
this  refractive  index  decrease  over  the  temperature  range, 
The  half-power  full-width  of  the  fundamental  transverse 
mode  dl  decreased  over  28-42  *C  in  a  manner  consistent 
with  the  Sij  calculated  for  the  second  mode.  The  second- 
mode  lobe-separation  angle  6^  was  chosen  for  calculations 
for  better  accuracy. 

An  interesting  feature  of  this  laser  is  the  cavity  selection 
of  modes  at  an  intermediate  temperature  of38  *C,  where  Fig. 
4  suggests  both  lower-order  modes  may  exist  dufing  a  cur¬ 
rent  pulse.  In  Figs.  8(d)  and  8(e)  the  far  field  is  filtered  with 
the  monochromator  and  the  output  of  the  photomultiplier  is 


fed  into  the  sampling  scope  with  a  50-/2  termination.  The  top 
pulse  of  each  picture  is  the  current  pulse  while  the  lower  is 
the  optical  pulse  of  the  laser.  It  was  established  during  the 
experiment,  from  output  pulses  as  in  Fig.  8  and  from  near¬ 
field  pictures,  that  the  laser  switched  modes  during  the  cur¬ 
rent  pulse.  The  discussion  to  follow  explains  this  multimode 
switching  with  current  magnitude  at  an  intermediate  tem¬ 
perature.  Theoretical  results  accompanying  the  experiment 


REFRACTIVE  INOEX  Vt 


HEAT  SINK  TEMPERATURE  (*C) 

FIG.  3.  Lasing  wavelength  vs  heat-sink  temperature  of  PL1 74- 1-31. 


HO.  7.  Experimental  (solid  tine)  vs  heat-sink  temperature  and  live- layer 

waveguide  calculated  6„  (dashed  line)  vs  active- region  index. 
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FIG.  8.  Time  display  of  input  laser  current  and  optical  output.  In  (a)  an 
input  current  pulse  is  drawn.  Jltli  and  Jtk:  are  the  threshold  current  densities 
of  the  first  two  modes,  respectively.  In  (b),  optical  flux  0,  at  the  fundamen¬ 
tal  mode  wavelength  is  shown  to  exist  in  regions  I  and  III.  In  (c),  the  second 
mode  optical  flux  02  is  detected  in  region  II  at  the  second  mode  wavelength. 
In  (d),  the  experimental  output  pulse  is  shown  in  the  lower  trace  at  the  fun¬ 
damental  mode  wavelength.  The  top  pulse  is  the  input  current.  In  (e),  the 
lower  trace  is  the  second-mode  output  it  its  wavelength. 


show  the  cavity  selection  is  sensitive  to  gain/mode  profile 
coupling  and  index  change  with  excess  carrier  injection. 

IV.  DISCUSSION 

The  data  from  Fig.  3  were  used  with  the  data  from  Ta¬ 
ble  I  for  a  multilayer  waveguide  analysis  of  the  transverse 
modes.  This  method  was  used  first  to  check  the  accuracy  of 
material  parameters.  In  Fig.  6,  the  dashed  line  shows  a  close 
fit  with  the  solid  measured  curve  of  the  far  field.  The  pip  on 
one  shoulder  is  a  leaky  wave  coupling  to  the  substrate  layer. 
As  previously  mentioned,  Snell's  law  dictates  a  fixed  refrac¬ 
tive  index  for  the  substrate.  It  was  noticed  in  the  experiment 
that  the  pip  remained  stationary  in  angular  displacement 
with  temperature  change  20-40  *C  and  with  current  drive. 
Also,  the  pip  remains  about  4  dB  below  the  main  lobes  over 
temperature  and  current  variations.  This  means  the  thick¬ 
ness  and  refractive  index  of  layer  3  in  Fig.  2  should  be  fixed  in 
order  to  hold  the  pip  stationary.  The  only  feature  to  change 
for  the  second  mode  is  the  angular  separation  of  the  major 
lobes  designated  as  These  data  are  represented  in  Fig.  7 
and  will  be  discussed  subsequently. 

The  only  way  to  match  the  0^  change  with  temperature 
is  to  slightly  decrease  the  active-layer  refractive  index  with 
increased  current  drive.  The  increase  in  threshold  current  of 
a  mode  with  temperature  occurs  due  to  a  modal  confinement 
change,  internal  quantum-efficiency  decrease,  or  a  carrier¬ 
lifetime  or  absorption-coefficient  change.20  Near-field  mea¬ 
surements  did  not  show  a  shift  of  the  transverse  mode  with 
temperature.  The  refractive-index  change  of  each  layer  with 
temperature Srf/STcan  be  shown  to  be  1.5x  10  */°K2'  For 
a  twenty-degree  range,  20-40  *C,  6i/~0.003.  On  the  other 
hand,  the  index  change  due  to  dispersion21  for  a  layer  is  8ij/ 

0.8/e V.  For  6A.  =  20  A  frdm  Fig.  5  (the  first  high-order 
mode),  it  can  be  shown  that  —  0.003/A. 

The  dispersion  and  heating  effects  on  refractive  index  of 
each  layer  tend  to  cancel.  Therefore,  with  the  evidence  that 
the  near  field  remains  constant  with  temperature,  the  pip  on 
the  shoulder  of  the  far  field  remains  stationary  and  disper¬ 
sion  and  heating  effects  tend  to  cancel,  it  is  assumed  that  the 
dominant  refractive-index  change  is  that  of  the  active  layer 


TABLE  I.  Typical  values  for  the  physical  constants  used  in  the  ambipolar 
diffusion  solution* 


Layer  I 

Layer  2 

Layer  3 

Doping 

p-HT'* 

P-Sxio11 

«-5  Xl0IT 

Al  concentration 

20% 

0 

2% 

V 

3.48 

variable 

3.5865 

Absorption 

5 

20.8 

2.4 

Width 

2.65x10-* 

.  .  . 

ft. 

3132 

3649 

3649 

138 

167 

167 

D. 

81.0 

198 

128 

Df 

3.67 

4.65 

4.33 

r 

7.69XI0-* 

3.59XIO-* 

8.46X10'* 

D* 

80.6 

156 

68 

L  * 

7.87x10-* 

2.37x10-* 

2.4x10-* 

*  All  quantities  are  in  cp  units,  with  the  exception  of  mobility  which  is  in  the 
mixed  units  cmVV  sec  and  absorption  which  is  in  cm" *. 
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TABLE  II.  Calculated  values  for  threshold  current  density,  modal  absorp¬ 
tion  and  internal  differential  quantum  efficiency  obtained  from  theambipo- 
lar  diffusion  solution 


Mode  1 

Mode  2 

Threshold  current  density  |kA/cm'i 

15.58 

1693 

Modal  absorption  (cm  \ 

19.74 

6.38 

Internal  differential  quantum  efficiency  foy  I 

66.8 

72.8 

due  to  injection  level.  The  data  of  Tables  I  and  II  and  Fig.  3 
are  used  to  find  the  solution  to  the  ambipolar  diffusion  equa¬ 
tion.  A  rise  in  temperature  is  represented  by  a  reduction  in 
refractive  index  as  required  in  matching  the  far-field  pattern 
of  Fig.  6. 

With  the  relevant  material  parameters  listed  in  Fig.  3,  a 
solution  was  found  for  the  ambipolar  diffusion  equation.  In 
Fig.  9  the  threshold  current  density  for  the  first  two  modes  is 
plotted  as  a  function  of  the  active-region  refractive  index. 
For  larger  i j2  values  (corresponding  to  a  larger  index  step) 
the  second  mode  has  the  lowest  threshold  due  to  superior 
coupling  to  the  excess  carrier  distribution.  But  as  the  active- 
region  index  is  lowered  and  the  second  mode  approaches 
cutoff,  its  coupling  efficiency  diminshes  and  the  fundamen¬ 
tal  mode  emerges  with  the  lowest  threshold. 

Typical  material  parameters  calculated  (fundamental 
mode  at  rj2  =  3.5948)  are  shown  in  Table  I.  Since  many  of 
the  parameters  depend  on  injection  levels,  the  quantities  in 
Table  I  vary  slightly  over  the  graphs  presented.  With  these 
calculated  values,  the  usefulness  of  the  ambipolar  equation 
can  be  shown  for  the  active  layer.  The  calculated  hole-injec¬ 
tion  efficiency  is  99.4%  and  the  electron-injection  efficiency 
is  86%  for  this  case.  This  calculation  predicts  a  crossover 
point  at  a  slightly  lower  tj2  value  than  the  multilayer  mode 
range  in  Fig.  7.  This  is  because  in  the  multilayer  model,  the 
presence  of  the  substrate  with  a  relatively  high  refractive 
index  increases  the  active-region  eigenvalue.  For  instance,  if 
the  buffer-layer  thickness  dx  were  decreased,  the  active-lay¬ 
er  eigenvalue  would  shift  even  more. 


FIG.  9.  Ambipolar  diffusion  solution.  Threshold  current  density  for  the 
flrst-and  second-order  modes  is  plotted  against  the  active-region  index  rj,. 
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Figure  9  is  useful  in  interpreting  the  phenomena  shown 
in  Fig.  4.  The  fundamental  mode  does  not  occur  below  30  *C. 
A 1 20  °C,  the  second  mode  has  a  threshold  current  density  of 
about  12  kA/cm2.  However,  abruptly  above  30  *C,  the  fun¬ 
damental  mode  appears.  From  Fig.  8,  a  slight  change  in  t/2 
on  the  order  of  Srj2  =  —  0.002  is  sufficient  to  cause  a  cross¬ 
over  in  threshold  currents  between  the  two  modes.  Both 
Figs.  4  and  9  show  that  the  second-mode  threshold  increases 
sharply  after  the  crossover.  Threshold  current  densities  pre¬ 
dicted  by  the  ambipolar  solution  are  in  fair  agreement  with 
experimental  values. 

Modal  output  power  as  a  function  of  current  density  is 
plotted  in  Fig.  10.  Because  the  linear  theory  of  gain  is  as¬ 
sumed,  only  one  mode  at  a  time  is  calculated.  At  power  levels 
above  threshold,  the  modal  gain  is  assumed  pinned  at  its 
threshold  value.  Note  that  the  second  mode,  despite  higher 
threshold  current  density  and  reduced  coupling  efficiency, 
has  a  higher  internal  differential  quantum  efficiency  than  the 
fundamental  mode.  Thus,  for  a  sufficient  overdrive  beyond 
threshold,  the  second  mode  will  eventually  be  the  mor;  effi¬ 
cient  overall  and  will  dominate.  An  explanation  for  this 
higher  internal  differential  quantum  efficiency  can  be  seen  in 
Table  II.  The  modal  absorption  of  the  second  mode  is  low 
because  a  significant  portion  of  its  energy  is  propagating  out¬ 
side  the  lossy  active  region  whereas  most  of  the  fundamental 
mode  is  well  confined  to  the  active  region.  Since  modal  loss 
dominates  slope  efficiency,  the  second  mode  has  the  higher 
value.  This  higher  value,  in  conjunction  with  a  higher 
threshold  due  to  poorer  coupling  efficiency,  causes  the  mode 
efficiency  crossover  point. 

With  these  results,  the  phenomena  of  mode  crossover  in 
this  laser  can  be  understood.  At  20  *C  the  second  mode  has 
the  lower  threshold  and  higher  efficiency  and  consequently 
only  the  second  mode  appears.  At  an  intermediate  tempera¬ 
ture,  the  second  mode  will  appear  at  high-enough  current 
overdrive  due  to  better  slope  efficiency.  Above  *2  *C,  the 
threshold  for  the  second  mode  becomes  prohibitive.  Figure  8 
depicts  a  time  resolution  between  the  two  modes  at  an  inter¬ 
mediate  temperature.  Figure  8(a)  represents  an  input  drive 
pulse.  The  two  threshold  current  densities  for  the  two  modes 
are  labeled  J,  and  J2.  As  the  current  pulse  rises  above/,  the 
first  mode  turns  on  in  region  I.  Output  flux  d>,  for  the  first 
mode  is  represented  in  Fig.  8(b).  As  the  input  current  pro¬ 
gresses  to  /2,  region  II  is  entered.  The  second  mode  turns  on, 
as  shown  in  Fig.  8(c),  as  <P2.  Finally,  with  the  drop  of  input 
current,  region  III  is  entered  and  first  mode  reappears.  As 
shown  in  Figs.  8(d)  and  8(e),  this  situation  is  observed  experi¬ 
mentally.  In  Fig.  8(d),  the  monochromator  is  set  on  the  first¬ 
mode  wavelength  and  the  lower  pulse  is  the  output  power. 
Reviewing  the  two  output  pulses,  it  can  be  seen  the  two 
modes  have  a  time  separation  during  a  current  pulse.  This  is 
explained  by  the  fact  that,  from  Fig.  9,  the  second-mode 
threshold  can  still  be  reached  slightly  after  crossover.  From 
Fig.  1 0,  the  second  mode  clearly  is  predicted  to  dominate  the 
first  mode  if  threshold  can  be  reached,  due  to  better  effi¬ 
ciency. 

Theoretical  mode  and  excess  carrier  coupling  for  the 
laser  is  shown  in  Fig.  1 1,  where  the  excess  carrier  profile  is 
plotted  across  the  active  layer  for  the  two  modes  operating  at 
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FIG.  10.  Am  bipolar  diffusion  solution. 
Output  power  it  plotted  against  threshold 
current  density  for  the  first  two  modes. 
The  slifiht  curvature  in  the  mode  I  curve  is 
caused  by  earner  redistribution.  The  in¬ 
dex  Vi  -  3  3948. 


CURRENT  DENSITY  J  (KA/CM*) 


two  power-output  levels  (threshold  and  0.67  w/mil.).  Includ¬ 
ed  on  the  carrier  plot  is  the  normalized  modal  distribution. 
Several  factors  should  be  noted.  First,  the  second-mode  dis¬ 
tribution  is  clearly  penetrating  the  passive  n  layer,  indicating 
a  near  cutoff  situation  and  a  modal  loss  dominated  by  the  n- 
layer  absorption.  With  all  curves  normalizd,  a  simple  pro¬ 
duct  of  Sp  and  E2,  integrated  over  the  active  region,  shows 
the  fundamental  mode  better  coupled  than  the  second,  as 
expected.  Second,  it  can  be  seen  that  the  carrier  profile  peaks 
on  the  side  where  holes  are  injected.  This  profile  contradicts 
the  notion  that  electron  injection  dominates  and  the  peak 
occurs  on  the  n  side.  Since  the  charge  cloud  is  ambipolar,  it 
must  accommodate  the  low-hole  mobility  to  maintain  quasi- 
charge  neutnJity.  Hence,  the  profile  is  skewed  toward  the 
hole  injection  or  p  side  of  the  active  layer.  As  expected,  car¬ 
rier  profiles  change  shape  with  current  density. 

As  a  final  note,  we  use  the  refractive-index  dependency 
of  the  active  region  on  temperature  and  hence  current  den¬ 
sity  to  estimate  the  free-carrier  effect.  The  influence  of  free 
carriers  on  the  refractive  index  in  a  material  is  determined 
from  the  classical  theory  of  dispersion.  The  refractive  index 
in  the  presence  of  free  carriers  rj  is  related  to  the  refractive 
index  in  the  absence  of  carriers  rj  by 


v-rfi.-Sfe). 


Here,  m*  is  the  effective  mass  of  the  carrier,  q  is  the  charge  of 
an  electron,  N  is  the  carrier  density,  a  =  2 irc/A,  and  c0  is  the 
vacuum  permittivity.  Since  the  change  in  index  is  much 
smaller  than  rj,  Eq.  (IS)  can  be  rearranged  to 

(161 

m*(o  rj 

In  other  words,  the  change  in  index  due  to  free  carriers  is 
directly  proportional  to  the  concentration  of  free  carriers 


KN,  (17) 

where  K  is  a  constant  we  wish  to  determine.  By  plotting 
refractive-index  change  versus  carrier  increase,  the  constant 
K  would  be  the  slope  at  a  particular  point 

d(6Vtc)  =  Kd(N)  (18) 

or 

(19) 

dm 

A  change  in  a  particular  parameter,  say  lifetime  r,  with  injec¬ 
tion  would  be  reflected  in  Eq.  (17)  having  a  nonzero  second 
derivative.  We  take  Srjfc  from  Fig.  7  and  conclude  that  this 
index  change  is  very  nearly  a  linear  function  of  temperature. 
However,  we  find  from  Fig.  4,  carrier  concentration  N,  relat¬ 
ed  to  the  peak  threshold  current,  is  not  quite  a  linear  function 
with  temperature.  Therefore,  under  the  assumption  the  free- 
carrier  effect  is  the  dominant  cause  of  observed  pattern 
change ,  the  constant  K  would  be  precisely 

R  d(&i,kydT 
d(N)/dT  ' 

This  expression  contains  the  information  that  AT  is  not  a  true 
'‘constant”  but  is  itself  a  slight  function  of  injection  level. 
Nevertheless,  it  is  instructive  to  estimate  K  using  Eq.  ( 19)  to 
compare  with  conventional  theory.  Using 

A r  = 

qdi 

and  a  straight  line  approximation  to  d  (AT )  over  the  tempera¬ 
ture  range  20-42  *C  we  get 

K  =  —  3.1  X  10-21  cmJ,  r  =  3  nsec 
and 

AT  =  —  6.2  X 10  ~  31  cm3,  r  =  2  nsec. 
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FIG.  1 1  Excess  carrier  densities  plotted 
across  the  active  region  for  the  ambipo- 
lar  diffusion  equation.  Iq  (a),  the  second¬ 
mode  distribution  is  plotted.  In  (b).  the 
first-mode  distribution  is  shown.  Note  in 
(a)  that  the  second  mode  “tail”  extends 
into  the  n  cladding. 
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These  results  can  be  compared  to  K  =  —  4X 10-21  cm3 
from  Thompson22  and  K  —  —  (5-1 1) X  10~21  cm3  from  Sel- 
way.23  Thompson  used  r=3  nsec  for  his  calculation.  Sever¬ 
al  other  values  for  AT  have  recently  been  reported.  Olsson  and 
Tang24  report  K  -  -  (4.9  ±  0.4)X  10-21  for  r  =  2.2  nsec. 
Henry  tt  al.13  report  K  =  —  (1.8  ±  0.4)  X 10" 21  for  lasers 
with  carrier  densities  estimated  at  1.02  X 10"  cm-3. 

V.  CONCLUSION 

A  solution  for  the  ambipolar  diffusion  was  presented 
for  the  active  region  of  a  double  heterojunction  laser.  Along 
with  the  solution  of  modes  from  Maxwell’s  equations  in  a 
multilayer  waveguide,  this  combined  description  was  ap¬ 
plied  to  a  broad-area  heterojunction  laser.  The  refractive- 


index  change  of  the  active  region  with  temperature  was  esti¬ 
mated  from  far-field  pattern  matching  of  the  transverse 
modes.  Transverse-mode-selection  conditions  were  predict¬ 
ed  with  the  diffusion  equation  and  compared  with  the  trans- 
verse-mode  selection  of  the  sample  laser. 

Mode  threshold  was  shown  to  depend  on  coupling  effi¬ 
ciency  between  excess  carriers  and  modal  distribution  in  the 
active  layer.  The  differential  internal  quantum  efficiency  de¬ 
pends  primarily  on  mode  absorption  and  injection  effi¬ 
ciency.  Mode  crossover  occurs  when  the  mode  with  the 
highest  threshold  also  has  the  highest  slope  efficiency. 

Carrier  transport  in  a  semiconductor  laser  is  an  ambi¬ 
polar  diffusion  process  which  can  be  described  using  the  re¬ 
laxation  time  approximation.  In  bulk  regions,  quasi-charge 
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neutrality  is  satisfied  because  of  the  strong  internal  field 
between  electrons  and  holes.  The  active  region  of  a  semicon¬ 
ductor  laser  undergoes  double  injection:  hole  injection  from 
the  p  side  and  electron  injection  from  the  n  side.  In  devices 
with  good  injection  efficiencies,  the  injected  carrier  profile  is 
higher  at  the  p-p,  rather  than  at  the  p-rt,  interface  because  of 
low  hole  mobility. 

The  device  modeled  here  shows  a  distinct  peak  of  carri¬ 
ers  on  the  p-p  interface.  The  slope  of  the  carrier  profile  in¬ 
creases  with  increased  current  density.  Since  these  carrier 
profiles  demonstrate  significant  skewing,  it  is  clear  that  this 
can  affect  mode  selection  if  the  modal  shape  couples  poorly 
with  the  carrier  concentration. 
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Abstract 

A  numerical  method  and  the  effective-index  method 
are  applied  to  a  three-layer,  constant  thickness  di¬ 
electric  waveguide  with  smoothly  varying  dielectric 
constant  Inside  the  active  layer  and  constant  permit¬ 
tivity  In  the  confining  layers.  The  results  of  the  two 
methods  are  compared  In  terms  of  the  propagation  cons¬ 
tant  >  calculated  by  each  method.  Application  of  the 
effective-index  method  facilitates  a  physical  under¬ 
standing  of  dielectric  waveguide  modes  as  well  as  pro¬ 
viding  an  efficient  approximate  method  for  calculating 
mode  behavior. 
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Introduction 

Several  papers  [1-3]  have  analyzed  mode  propagation 
In  dielectric  waveguides  with  spatially  varying  refrac¬ 
tive  Index,  usually  approximated  by  parabolic  or  tanh^ 
functions,  which  go  to  infinity  at  large  distances  from 
the  reference  point  x  ■  0.  We  will  'iae  the  approxima¬ 
tion  [4] 

X  "  *s  +  Ax/cosh2  (x/Xq)  (1) 

to  describe  the  variation  of  < .  This  la  In  closer  cor¬ 
respondence  with  the  physical  situation  and  leads  to 
equations  with  known  solutions.  The  disadvantage  In 
this  case  Is  the  fact  that  the  field  solutions  consist 
of  a  finite  (possibly  empty)  set  of  confined  (trapped) 
modes,  an  Infinite  set  of  diverging  "leaky"  modes  and 
a  continuum  of  solutions  that  will  be  designated  as 
"radiation"  modes. 

The  2-dlmenslonal  character  of  the  problem  due  to 
the  variation  of  the  refractive  Index  In  the  lateral 
(x)  and  transverse  (y)  directions  requires  the  uee  of 
approximate  methods  or  numerical  solutions.  Among  the 
former,  the  effective  Index  method  Is  the  most  popular 
([2],  [4]  -  [8]).  Direct  numerical  Integration  Is  pos¬ 
sible,  but  the  computation  time  required  is  much  larger 
than  that  needed  for  other  methods  like  the  one  used 
In  [1],  which  we  will  apply  in  this  paper. 

The  following  sections  describe  the  class  of  wave- 
guldes  considered,  the  numerical  method  and  the  appli¬ 
cation  of  the  effective- Index  method  to  our  problem, 
ending  with  conclusions. 


Description  of  Structures 


Figure  1  shows  the  structure  considered  In  this  pa¬ 
per.  The  confining  layers  A  and  C  are  assumed  identi¬ 
cal,  their  refractive  indices  being  described  by 
2 

'a  ■  «C  ;  ■  1  QAnA/k0  (2> 

describes  the  power  loss  In  these  layers  and  Is 

constant  with  distance.  The  active  layer  has  a  constant 
thickness  d,  with  a  refractive  Index 

*(x)  •  n2(x)  -  1  a(x)  n(x)/kn  (3) 


whose  dependence  with  x  la  conaldered  to  be  reaaonably 
well-approximated  by  (1).  In  that  equation. 


*0  ■  “s 

(4a) 

x(0)  -  nj  +  1  gQ  n0/kQ 

(4b) 

"s  "  1  as  ns/ko 

(4c) 

Xq  Is  a  parameter  related  to  the  width  of  the  stripe 

contact  In  the  case  of  a  semiconductor  laser.  The 
values  of  power  attenuation  coefficient  and  refractive 
index  Inside  the  active  region  far  away  from  the  stripe 
are  c»s  and  ng,  respectively.  The  quantity  gQ  repre¬ 
sents  the  peak  power  gain  under  the  stripe  (x  -  0) , 
where  the  refractive  Index  Is  n^.  For  An  «  nQ-ns>0 

the  mode  will  be  lndex-gulded,  while  for  An  <  0  It  will 
be  Index  antl-guided,  and  this  latter  effect  can  even¬ 
tually  offset  the  guiding  effect  of  the  gain  distribu¬ 
tion. 


Numerical  Solution 

We  follow  the  method  used  In  [1].  Maxwell's  equa¬ 
tions  are  solved  for  both  the  active  and  the  confining 
layers,  requiring  as  boundary  conditions  that  the 
general  solutions  Inside  and  outside  the  active  layer 
and  their  normal  derivatives  match  at  the  Interfaces 
y  *  +  d/2.  We  also  demand  chat  these  solutions  vanish 
at  x  »  +  «°,  y  •  +  •>.  Application  of  the  method  of 
separation  of  variables  results  In  a  vertical  field 
solution  Inside  the  active  region  of  the  form 


*(y)  - 


and  a  differential  equation 

+  ^k0  ics+’,2~q2+k0 

dx 


coah  (x/xQ) 


where  q  Is  the  separation  constant.  Equation  (6) 
possesses  solutions  of  the  form 

t-b0  b^-l-dj 

*t(x)  *  [cosh(x/x0)l  (tanh[x/Xgl ) 


where  C*(«)  *ra  Cegenbauer  polynomial#  (9)  and 
b0  -  -  *5+  &  ♦  k2  Xq  Ak)^ 

To  sstlsfy  the  boundary  conditions  we  require 


Re{b0  -  t)  >  0 


For  Refbg)  >  0,  the  fundamental  mode  (l“0)  Is  trapped. 
Modes  of  order  1  such  that  Relb^  -  t)  <  0  are  still 

solutions  of  (6)  but  diverge  as  jx|  -  “  and  will  be 
designated  as  "leaky".  Radiation  modes  would  be  des¬ 
cribed  by  other  solutions  of  (6)  for  srbltrary  (non- 
integer)  eigenvalues  [9].  The  general  solution  will 
be  s  linear  combination  of  the  few  discrete  trapped 
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inodes  (8)  satisfying  (9),  plus  an  Integral  over  the 
continuum,  and  the  leaky  modes  must  be  excluded  If  the 
field  has  to  vanish  for  |x|  -»  *>.  The  Importance  of  the 
continuum  can  be  expected  to  decrease  as  the  number  of 
trapped  modes  Increases.  Since  this  number  is  rela¬ 
tively  large  for  structures  with  modal  gains  that  are 
high  and  not  very  sensitive  to  the  refractive  Index 
step  An  *  Oq  -  n,.,  this  continuum  will  be  neglected. 

Hence,  the  general  even  solution  Inside  the  active  re¬ 
gion  Is  approximated  by 

£ 

max 

Vx*y)  s  I  A«  c°s(qt  y) 

*-0,2....  1  (10) 


where  £  is  the  maximum  even  value  of  £  for  which 
max 

l£(x)  is  confined.  Solution  of  the  wave  equation  for 

the  confining  layers  proceeds  as  In  (1).  Matching  the 
functions  and  their  derivatives  with  respect  to  y  at 
the  boundary  y  »  d/2  and  applying  the  orthogonality 
relation  for  the  trapped  modes  results  in  a  finite  sys¬ 
tem  of  homogeneous  linear  equations 

(ST  -  I)A  -  0  (11) 

-T 

where  A  =  (Ag,  k^,  ....  A—  ) ,  I  is  the  unit  matrix, 

max 

and  the  matrix  elements  of  0  are  given  by 


4  cos(q£  d/2)  *'n^ 
rx0  q£,  sin(qtl|)  >'s~, 


f  _  —  2  2  2  k 

where  I£  £,(y)  *  j  (x)  (x)  (x  ”1  _ka)  d* 

(13) 

!lif(>.)  being  the  normalized  Fourier  cosine  transform  of 
tl£(x),  and  £,£'  «  0,  2,  ....  Imj|.  The  q£  satisfy 

"  V2  +  kg  xg  +  (bQ  -  1)2/Xg  (14) 


a  parameter.  We  notice  that,  for  each  value  of  gQ,  An 

can  be  decreased  up  to  a  certain  value  beyond  which  the 
fundamental  mode  becomes  leaky. 

Effective-Index  Calculation 

The  effective-index  method  consists  basically  of 
reducing  a  two-dimensional  problem  to  an  equivalent 
one-dimensional  one.  In  our  case,  the  two-dimensional 
character  of  the  problem  is  given  by  the  dependence  of 
the  dielectric  constant  on  x  and  y.  As  a  first  appro¬ 
ximation,  the  variation  in  one  direction  (in  our  case, 
x)  is  neglected;  this  la  justified  if  this  variation 
is  much  less  than  that  in  the  y  direction.  This  is 
equivalent  to  approximating  the  waveguide  with  a  sim¬ 
ple  3-layer  guide  whose  dielectric  constants  do  not 
vary  with  x.  The  solution  of  this  problem  yields  the 
transverse  or  vertical  variation  of  the  field.  Next, 
the  original  equation  describing  the  2-dlmenslonal 
equation  in  x  can  be  solved  for  the  lateral  variation 
of  the  field,  and  the  overall  solution  is  approximated 
by  the  product  of  this  lateral  field  and  the  vertical 
field  found  from  the  3-layer  problem.  We  start  with 
the  wave  equation  in  2  dimensions: 


V2  f  +  [y2  +  k2  x (x,y) ]  *  -  0 

(15) 

For  the  simple  3-layer  guide  we 

assume  x(x)  -  >. 

inside 

the  active  layer.  Now,  we  transform  (15)  making 
t(x,y)  -  it' (x)  »(y),  multiplying  it  by  **(y)  and 
grating  it  over  y  from  -«>  to  ■»,  and  obtain 

inte- 

liv  +  ,2  2  +k2  K  +k2 

Ax  .. 

_  f  _  j 

/I  £.  \ 

.  2  lY  qeff  k0  S  eff  k0 

dx 

2 

cosh  (x/Xg) 

With  q2ff  -  r  k2  x0  -  p2  - 

k0  *A 

(17a) 

*S  eff  *  F  ‘s 

(17b) 

Ax  »  r  Ax 

eff 

(17c) 

The  system  of  equations  (11)  has  a  non-trivial  solution 
only  if  det (T.-I )  =0 .  Numerical  computation  of  the  roots 
yields  the  possible  values  of  the  propagation  constant  y. 

Numerical  Results 

The  method  was  applied  to  a  structure  described  by 
the  following  parameters: 

nA  -  3.38,  aA  •  50  cm  1,  •  50  cm  1 ,  gg  -  200cm  * 

ng  •  3.595,  Xg  -  6pm. 

Direct  solution  of  (11)  involves  computation  of  the 
n{8,.  In  general,  £1^,  -  but  1££,  -  T£,r  Values 

of  the  modal  loss  *  Re(y)  were  computed  by  evaluating 
I £ £ ,  (> )  using  the  exact  expression  (13).  This  has  to 

be  done  once  for  every  value  of  y.  Instead  of  solving 
the  complete  system  (11),  we  Btart  with  a  lxl  matrix, 
go  on  to  a  2x2  matrix,  etc.,  and  observe  that  the  re¬ 
sult  converges  relatively  fast  for  a  2x2  case,  which 
is  considered  suffic)  'nt  for  our  approximation.  In 
spite  of  this,  the  computation  time  is  still  impracti¬ 
cal  ly  high.  An  increase  in  speed  by  nearly  two  orders 
of  magnitude  while  still  maintaining  good  accuracy  is 
achieved  expanding  the  radical  in  (13)  using  the  bino¬ 
mial  theorem  (which  results  in  integrals  that  do  not 
depend  on  >)  and  retaining  the  first  three  terms. 

Figure  2  shows  plots  of  the  modal  loss  »  Re(y)  as  a 
function  of  An  with  gg,  the  gain  under  the  stripe,  as 


where  T  is  the  filling  or  confinement  factor. 

Equation  (16)  has  the  same  form  as  eq .  (6).  It  will 
also  have  polynomial  solutions  similar  to  (  7  )  that 
represent  trapped  modes: 

£-b  eff  b  eff-l+k 

*j(x)  -  [cosh(x/xg)l  Ct  (tanh  — ) 

>.18 ) 

£  -  0,  1,  2,  ... 

where  bg  ^  is  defined  as  in  (8)  with  Ax  replaced  by 

F  Ax  and  p  is  the  eigenvalue  of  the  simple  3-lavcr  problem 
For  the  fundamental  mode,  £  -  0,  and  we  obtain 

y2  ■  *k0  ‘A  -  (b0  eff/x0>  -  (”2-k0  r  AO  (19) 
for  the  propagation  constant. 

Discussion 

Values  of  y  obtained  using  (19)  are  compared  with 
those  obtained  with  the  numerical  method  in  fig.  3. 
Solutions  are  very  close  for  all  values  of  An  for  which 
the  mode  exhibits  a  gain  which  is  relatively  high  ad 
with  low  sensitivity  to  An.  The  results  differ  most  in 
the  range  of  An  for  which  the  mode  has  a  net  loss  or 
has  a  relatively  low  gain  with  higher  sensitivity  to  An. 

Figure  4  gives  the  required  value  of  the  peak  power 
gain  gg  under  the  stripe  to  obtain  a  given  modal  gain 

C.  as  a  function  of  An,  using  the  effective-index  method. 
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Also  Included  Is  the  region  for  which  the  fundamental 
mode  becomes  leaky.  The  vertical  confinement  factor  T 
did  not  vary  appreciably  with  An;  a  typical  value  for 
the  case  considered  was  r  =  0.4963.  We  see  that 
lateral  variations  in  the  refractive  index  affect  the 
gain  much  more  by  altering  the  lateral  field  distribu¬ 
tion  than  by  affecting  the  vertical  variation. 

The  effective-index  method  is  seen  to  be  a  fast  and 
relatively  accurate  way  to  obtain  the  field  distribu¬ 
tions  for  Lhe  class  of  waveguides  considered,  for  which 
the  numerical  method  we  used  is  not  practical  for  ex¬ 
tensive  modeling  due  to  its  long  computation  time  in 
spite  of  all  approximations,  The  remarkable  agreement 
between  the  ef fective- index  method  and  experimental 
results  found  by  other  workers  ((10])  increases  our 
confidence  in  this  powerful  approximate  method. 


confining  loyw 
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Figures 

1.  Waveguide  structure  considered 

2.  Modal  loss  as  a  function  of  An  with  peak  power  gain 
gfl  as  a  parameter.  Regions  at  left  of  vertical 

lines  correspond  to  leakv  fundamental  mode. 

3.  Modal  loss  Re'v)  as  a  function  of  An  for  nuneric 
(lxl  matrix)  and  ef fective- index  methods. 

4.  Peak  power  gain  under  contact  stripe  as  a 

function  of  An  with  modal  gain  as  a  parameter. 
Shaded  area  shows  region  corresponding  to  leakv 
fundamental  mode. 


Figure  4 
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Modal  Solutions  of  Active  Dielectric  Waveguides 

by  Approximate  Methods 

A.  LINZ,  MFMBf'R.  1EEE,  AND  J.  K.  BUTLER,  SENIOR  MEMBER 


’  Abstract — Approximate  mvlhodv  are  uvH  lo  ohfain  the  rwvM  gropertlex 
of  xtripe- contact  semiconductor  Injection  lasers  using  a  planar  thiee -laser 

•  waveguide  model.  The  central  active  laser  has  a  dielectric  constant  that 
i  aries  smoothly  along  the  direction  parallel  to  the  helerojuncrion  boundaries. 
I  he  c, implex  dielectric  constant  under  the  stripe  contact  Is  dependent  on 
the  gain  and  approaches  a  constant  value  at  large  lateral  distances.  The  two 

.  methods  are  compared  In  terms  of  their  modal  propagation  constants.  An 
application  of  the  effective  Index  method  facilitates  a  physical  understand 
■  ing  of  dielectric  waveguide  modes  as  well  as  providing  an  efficient  calcttla- 

•  lion  procedure. 

I.  Introduction 

NALYSIS  OF  mode  propagation  in  dielectric  wave¬ 
guides  with  a  spatially  varying  refractive  index  Ins 
been  the  subject  of  several  papers  [  1 1-(3J.  Typically,  the 
variation  of  the  dielectric  constant  with  distance  has  been 
approximated  with  a  parabolic  profile  (1),  [2J  or  a  function 
of  the  form  tc  =  -  icn  +  ft3  tanh?( jc/jc0 )  |3).  Both  approxi¬ 
mations  have  the  disadvantage  that  the  value  of  k  goes  to 
infinity  at  large  distances  from  the  point  x  -  0,  which 
corresponds  to  the  axis  of  lateral  symmetry  of  the  struc¬ 
ture.  In  the  case  of  a  semiconductor  laser,  this  corresponds 
lo  the  region  below  the  center  of  the  contact  stripe.  Another 
approximation  that  eliminates  this  disadvantage  is  the  use 
of  a  function  of  the  form  (4) 

k  -  ks  +  Ait/cosh7(jc/.x0)  (1) 

lo  describe  ihe  variation  of  tc.  This  is  in  closer  correspon¬ 
dence  with  the  phvsical  situation,  since  tt  now  acquires  the 
value  for  r  :*•  ,t0.  Even  if  this  particular  form  of  vari¬ 
ation  of  k  does  not  describe  the  actual  variation  very 
closely,  it  retains  the  most  important  features,  and  leads  to 
equations  with  known  solutions.  The  disadvantage  in  this 
case  is  the  fact  that  the  field  solutions  consist  of  a  finite 
possibly  empty)  set  of  confined  trapped  modes,  an  infinite 
i 
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set  of  discrete,  diverging  "leaky”  modes,  and  a  continuum 
of  solutions  that  will  be  designated  as  "radiation”  modes, 
as  opposed  to  an  infinite  set  of  discrete  trapped  modes 
only,  as  in  the  parabolic  and  tanhJ(x/x0)  profiles.  Mode 
analysis  is  a  two-dimensional  problem,  since  the  refractive 
index  varies  in  both  lateral  (x)  and  transverse  (y)  direc¬ 
tions.  Therefore,  numerical  or  approximate  methods  need 
to  be  applied.  The  most  popular  and  effective  approxima¬ 
tion  method  is  the  "effective-index”  solution,  whereby  the 
two-dimensional  problem  is  reduced  to  an  equivalent,  one¬ 
dimensional  one  [2],  14J-|8J.  Numerical  methods  have  also 
been  developed.  For  example,  in  [1]  the  parabolic  variation 
is  used.  Maxwell’s  equations  are  solved  both  for  the  active 
layer  and  the  confining  layers,  and  then  superposition  is 
applied  to  both  types  of  solutions  to  form  a  general  expres¬ 
sion  for  the  field.  These  solutions  and  their  derivatives  are 
matched  at  the  boundaries  of  the  active  layer,  yielding  an 
infinite  system  of  linear  homogeneous  equations,  whose 
solutions,  numerically  obtained,  are  the  expansion  coeffi¬ 
cients  for  the  mode  in  terms  of  the  eigenfunctions  of  the 
active  layer  problem.  Of  course,  direct  numerical  integra¬ 
tion  of  the  two-dimensional  wave  equation  is  possible,  but 
the  computation  times  are  long  compared  to  those  required 
by  the  algorithm  discussed  in  this  paper. 

For  the  type  of  variation  considered  here,  a  general  field 
in  the  active  layer  must  be  expressed  as  a  superposition  of 
the  few  confined  discrete  modes  plus  an  integral  over  the 
continuum.  Leaky  modes  cannot  be  included  in  the  expan¬ 
sion  if  the  field  is  to  decrease  to  zero  for  large  distances 
from  the  stripe. 

Direct  application  of  the  numerical  method  used  in  (I] 
results  in  a  finite  set  of  linear  equations  (due  to  the  finite 
number  of  trapped  modes)  coupled  with  an  integral  equa¬ 
tion  (due  to  integral  over  the  continuum).  For  the  case  in 
which  only  one  trapped  mode  exists  (the  fundamental 
mode),  an  integral  equation  results,  which  can  in  principle 
be  solved.  However,  these  cases  will  be  seen  to  correspond 
to  structures  with  net  modal  loss  or  low  gains  very  sensitive 
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Leaky  modes  (/>  /max)  cannot  be  used  when  expressing  a 


general  confined  field  as  a  superposition  of  modes.  This 
results  in  a  finite  discrete  eigenvalue  spectrum,  and  the 
need  arises  for  a  continuum  in  order  to  have  a  complete 
orthogonal  set  of  functions.  We  will  seek  an  approximate 
solution  neglecting  the  continuum.  For  the  trapped  modes, 
the  product  solutions  will  be  of  the  form 

t,(x)  cos  (q,y) 

and  the  general  solution  will  be  approximated  by 

%(*,y)s  £  '4/'M*)cos(<7,>’) 

1-0.2.  ■  , 

(24) 

(25) 

for  a  field  even  in  x ,  where  /mai  is  the  maximum  even  value 
of  /  for  which  ipi  is  confined.  The  trapped  modes  satisfy  the 
orthogonality  relation 

7-  f°°  'h(x)'t'i(x)dx  =  N,Sl 

•"‘o' -00 

(26) 

The  first  few  normalization  constants  are 

,, 

°  r(;.0  +  !) 

(27a) 

_  2(A„-r)’-'3  r(M 

1  (*0-1)  r(b0  +  i) 

(27b) 

v  2yfn(b0  —  3/2)2( b0  —  \  )T(b0) 

2  (ft0-2)r(*o  +  i) 

(27c) 

where  F(  z )  is  the  factorial  function.  An  analytic  but  lengthy 
expression  for  the  TV,  exists  but  is  not  given  here. 

Solution  of  (R)  for  the  confining  layers  proceeds  as  in  [1 J. 
Separation  of  variables  is  applied  to  (8);  the  lateral  solu¬ 
tion  is  of  the  form  ccs(x.r)  whereas  the  vertical  solution  is 
a  decreasing  exponential.  Superposition  is  then  applied; 
this  results  in  an  integral  since  no  boundary  conditions  are 
available  that  would  result  in  a  discrete  spectrum.  The 
general  solution  will  be 

%(*')’)**/'  B(x)cos(x*) 

■'0 

exp[(x2  -  y2  -  kl)W\d/l- y)]  d\  (28) 

where 

k1  -  kZK 

Ka  h0KA- 

Now,  (25)  and  (28)  and  their  normal  derivatives  have  to 
be  matched  at  y-±d/ 2,  the  boundaries  of  the  active 
layer.  Matching  tpa  and  results  in  an  expression  for  B(x) 

•y  Lit  /  j  \ 

fl(x)=-  £  ^,/tycos  U-rU/tx)  (29a) 

t  •  0.2,-  •  .  V  ' 

where 

6(x)  =  /  ^4~  cos(Xx)dx.  (29b) 

J0  y'A', 


Matching  the  derivatives  respect  to  y  at  y  =  d/2  any 
applying  the  orthogonality  relation  (26)  yields,  after  using 
(29a) 

£  A^cosi  qd\lu{y) 

WJC0/_o,2. ....  '  21 

=  /f,.?,.sin(<7,.|)y/^  (30) 

with  /,  /'=  0,2,-  •  ,/m„  and 

//./(y)  =  /  'F/(x)4’/(x)(x2-72-^fl)l/2rfX  (31) 
Jo 

ilx), _L  r  (-QT(».fb5)... 

'  /A4..0  r(6„-l  +  i)m!(l-2m)! 

r’f(-i)*(i-")i  1 

•  £  {7-7 - r - t{x,b0+lk-l,xQy, 

(/even)  (32) 


where 


<(x.  v  «„)  -  Ti^r(x + i2r)r(x -,3?)  (?s| 

This  is  a  finite  system  of  linear  homogeneous  equations  of 
the  form 

(nr-/)yf=  0  (34) 

where  AT={A0.Ay,---,Af^),  I  is  the  unit  matrix,  and 
the  matrix  elements  of  SI  are  given  by 

Acos{q,d/2){N, 

My)  - 7-77 — My). 

irx0<7,'Sin(<7,.-  j/AV 
The  q,  satisfy,  from  (12a)  and  (12b) 


(35) 


g/2  =  Y2  +  *o*s  +  —  J- 


(36) 


The  system  of  equations  (34)  has  a  nontrivial  solution  only 
if  del(fl-/)  =  0.  Numerical  computation  of  the  roots 
yields  the  possible  values  of  the  propagation  constant  y. 

IV.  Numerical  Results 

The  method  was  applied  to  a  structure  described  by  the 
following  parameters: 

nA  =  3.38  g0  =  200  cm" 1 

0^  =  50  cm"’  n0  =  3.595 
a^SOcm"1  x0  <=  6  pm. 


Direct  solution  of  (34)  involves  computation  of  the  In 
general,  fi„.  *  R, ,,  but 

Values  of  the  modal  loss  =  Re(y)  were  computed  first  b> 
evaluating  /„.( y)  using  the  exact  expression  (31).  This 
required  one  computation  of  this  integral  for  every  value  of 
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boundaries  y  =  ±  d/2,  and  demand  that  and  de¬ 
crease  to  zero  as  x  -»  ±  oo,  and  also  that  -* 0  as  |/|  -» oo. 
These  boundary  conditions  are  homogeneous,  so  we  can 
apply  separation  of  variables. 

First  we  consider  the  solution  inside  the  active  layer 
making 

Substitution  in  (7)  results  in  a  vertical  solution  of  the  form 


-WI) 


"here,  from  (12a) 


9  |  V 

*  tt 


s  * 

CC  I  : 

O  - 1 — V 


and  a  differential  equation 

— +  L2k  +  v2  a 2  I  k1  lit  -  0 

dX 2  f  *°“'  +  Y  9  °cosh2(xAo)  *~° 


where  q  is  the  separation  constant.  The  substitutions 

(=tanh(xAo)  (11a) 

*-(l-<»)*/2IF(«)  (lib) 

M6o  +  1)  =  *o*oak  (12a) 

B2  =  {q2  ~y*  -  klics)x'o  (12b) 

in  (10b)  result  in 

(l-€J)Hw-(2X  +  l)|ir+a(a+2X)fF-0  (13) 

with 

X  =  R  +  )  (14a) 

«  =  b0-X  +  ).  (14b) 

Several  solutions  exist  for  (13).  When  a  is  an  integer 

/  =  0, 1,2,-  •  •,  we  have  polynomial  solutions,  which  are  the 
ultrasphcrical  or  Gegenbauer  polynomials  Cx(£)-  Physi¬ 
cally,  they  give  rise  to  an  infinite  number  of  discrete 
modes,  of  which  some  may  be  trapped.  The  radiation 
modes  would  be  described  by  other  solutions  that  satisfy 
(13)  for  arbitrary  a  (see  (9]). 

We  will  now  examine  the  polynomial  solutions  JF(£)  = 
C,x(()  in  detail.  For  a  an  integer  / 

X  =  60-/  +  )  (15a) 

B,~b0-I,  /  =  0, 1,2,  •  •  •  .  (15b) 

The  C,x(£)  are  defined  by 

C0X(£)  =  1  (16a) 

C,x(£)  =  2X£  (16b) 

(/+1)C,X  ,(()  =  2(H  X)^;(0-(/  +  2X-l)Cx_,(0- 

(16c) 

’'sing  (11a)  and  (11b)  we  obtain  the  expression  for  the 
modes 

fA*)”  [<™h  (  X  Ao )] '  ’  h°C Ihn  - ' + 5  ( tanh  ( x  A0  J ) 


DIELECTRIC  STEP  RATIO  A/A* 

4  Confinement  curves  (<&  versus  A'/A")  for  several  mode  orders. 
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*o*=-i  +  (i  +  AoJroAK)'/2’ 


3  12  3  4  5 

LATERAL  DISPLACEMENT  X/Xo 

Fig.  5.  Field  amplitude  versus  lateral  displacement  for  modes  of  order  0 
and  I  for  the  3  cases  shown  in  Fig.  4.  The  curves  for  /  -  I  are  magnified 
five  times. 


For  x  :*»  x0,  (17)  becomes 

tf'z(x)  “  1/2  exp  [_(xAo)(fy>~  0],C*0-,+  ,/2(l)- 

(19) 

So  to  satisfy  the  boundary  conditions,  we  require 

Re(&0-/)>0.  (20) 

For  a  given  b0,  if  Re(fc0)  >  0  there  exists  at  least  one  mode 
that  decreases  as  x  -♦  ±  oo,  for  /  =  0  (the  fundamental 
mode).  Modes  of  order  /  such  that  Re(fr0  -  /)  <  0  are  still 
solutions  of  (10b)  but  diverge  as  |x|  -*  oo,  and  will  be 
designated  as  “leaky”.  Notice  that,  if  R c{b0)  <  0,  even  the 
fundamental  mode  becomes  leaky.  Condition  (20)  can  be 
expressed  in  terms  of  the  real  and  imaginary  parts  of  Ak 
defining  the  quantities 

8  =  A'/A"  (21a) 

4>  =  2*0x0(A"),/2  (21b) 

where  Ak  =  A'+  /  A".  We  obtain 

m+i)-r  (22) 

4>2(2/  +  l) 

as  the  region  for  lateral  confinement.  Fig.  4  shows  4>  as  a 
function  of  8  for  several  values  of  /.  Fig.  5  shows  the 
behavior  of  the  |^|  for  sets  of  parameters  resulting  in  the 
(S,  <P)  points  marked  in  Fig.  4.  In  particular,  notice 
the  behavior  of  ,|  when  we  go  from  “1”  to  “2”  to  “3”. 
The  <P/(X)  diverge  for  />  lmax,  where 

/m„  =  /Nr[Re{b0)]. 


(18) 


(23) 
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to  the  dielectric  step  size.  For  structures  that  exhibit  higher, 
stable  gain,  several  discrete  trapped  modes  exist,  and  we 
achieve  reasonable  convergence  with  the  first  few  modes, 
so  that  considering  the  continuum  is  not  necessary.  We 
also  apply  the  effective-index  method  and  compare  it  with 
the  approximate  numerical  method  in  terms  of  the  prop* 
gation  constant  y,  which  is  calculated  as  a  function  of  the 
dielectric  step  size  in  the  active  layer  An  =  n0  -  ns,  where 
n0  and  ns  are  the  values  of  refractive  index  under  the  stripe 
and  far  away  from  it.  The  results  of  the  two  methods 
practically  coincide  for  the  cases  in  which  several  trapped 
modes  exist.  They  differ  appreciably  only  for  that  range  of 
An  for  which  only  one  trapped  mode  exists.  The  dis¬ 
crepancy  may  be  possibly  due  to  the  continuum,  but  this 
paper  does  not  investigate  this  matter  further.  The  follow¬ 
ing  sections  will  consist  of  a  description  of  the  class  of 
waveguides  considered,  followed  by  a  description  of  the 
approximate  numerical  method,  ending  with  the  applica¬ 
tion  of  the  effective-index  method  to  this  problem  and 
conclusions. 


II.  Description  of  Structures 

Fig.  1  shows  the  structure  considered  in  this  paper.  The 
confining  layers  A  and  C  are  assumed  identical,  their 
refractive  indices  being  described  by 

KA  =  Kcsnl-iaAnA/k0.  (2) 

aA  describes  the  power  loss  in  these  layers,  and  is  constant 
with  distance.  The  active  layer  has  a  constant  thickness  d, 
with  a  refractive  index 

K(x)  =  n2(x)-ia(x)n(x)//c0  (3) 

whose  dependence  with  x  is  considered  to  be  reasonably 
well  approximated  by 

ic(x)  =  ks  +  A(c/cosh2(x/x0)  (4) 

where 

Ak  =  k0  -  ks  (5a) 

*o  =  K(°)  =  "o  +  igono/k0  (5b) 

Ks  =  "s-,£l[s"s/*o-  (5c) 

x0  is  a  parameter  related  to  the  width  of  the  stripe.  The 
values  of  power  attenuation  coefficient  and  refractive  index 
inside  the  active  region  far  away  from  the  stripe  are  as  and 
r?s,  respectively.  The  quantity  g0  represents  power  gain 
under  the  stripe,  where  the  refractive  index  is  n0,  and 
k0  =  2v/\. 

Using  (4)  and  (5),  we  can  obtain  expressions  for  the 
variation  of  the  refractive  index  n  and  the  loss  a  (or  gain 
-a)  as  a  function  of  distance.  Figs.  2  and  3  show  n(x) 
versus  x/xn  and  o(jc)  versus  x/.t0.  For  A  n  =  n0  -  ns>  0. 
the  mode  will  be  index-guided,  while  for  An  <  0,  it  will  be 
index-antiguided.  In  this  latter  case,  it  will  still  be  confined 
because  of  the  gain  distribution,  but  the  field  will  be  more 
spread  and  the  m<">dnl  pain  will  be  low  (eventually  we  may 
have  a  net  power  lo  st  If  A n  is  negative  enough,  the 
guiding  effect  is  lost  and  the  modes  become  leaky.  The 
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Fig.  1.  Waveguide  structure  considered. 
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Fig  2.  Refractive  index  of  active  layer  versus  lateral  displacement  for 
as  -  50  cm~ ',  -  200  cm"  n0  —  3.595  for  different  values  of  ns. 


LATERAL  DISPLACEMENT  X/Xo 


Fig.  3.  Power  gain  -  o  for  active  region  versus  lateral  displacement  for 
as  -  50  cm" ',  n0  -  3.595,  ns  -  3.6  for  different  values  of  g0  ,-i- 

ix. 


condition  Art  >  0  results  in  strong  confinement  and  high 
and  stable  values  of  modal  gain. 

III.  Numerical  Solution 
We  assume  an  electric  field  of  the  form 

£*”*(*.>')•  (6) 

Following  [1],  we  apply  Maxwell’s  equations  to  the  struc¬ 
ture  in  Fig.  1,  and  obtain  inside  the  active  layer 

V(I'kfr  +  [yz  +  ^K(x)]Skfr  =  0  V) 

and  outside  the  active  layer 

V(j4,o  +  [y:  +  ^K4]'I'(,  =  0  (8* 

where  x(.r)  is  given  by  (4).  We  require  the  functions 

and  their  normal  derivatives  to  be  continuous  at  the 
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y,  i.e.,  every  iteration  in  the  solution  of  (34).  This  resulted 
in  unacceptably  high  computation  times,  in  part  because 
involves  repeated  use  of  a  gamma-function  routine  and 
also  because  the  integrand  is  an  oscillating  function  of  x- 
In  order  to  improve  speed,  instead  of  solving  the  complete 
system  (34)  we  start  with  a  1  X  1  matrix,  go  on  to  a  2x2 
matrix,  etc.,  and  observe  the  convergence  of  the  result  Fig. 
6  shows  Re(y)  versus  An  for  1  x  1  and  2x2  matrices.  \V> 
see  that  the  result  converges  relatively  fast,  so  a  2x2 
matrix  is  sufficient  for  our  approximation.  This  holds  for 
the  case  >  2.  For  lmax  =  0,  we  are  limited  to  the  lxl 

case.  In  Fig.  6,  for  An  < - 0.01, /m„  =  0.  For  An  > 

-0.01,  fmtx  increases  as  shown  in  the  figure. 

In  spite  of  this,  even  for  fmax  =  0,  the  computation  time  is 
impractically  high.  An  increase  in  speed  is  achieved  recog¬ 
nizing  that  (31)  can  be  written  as 

'Wy)-(-t,-*2)'/T 

t  ^  -  OO 

I  • 'F/(x)'fr(x)dX •  (37) 

For  common  values  of  k\  and  y,  |xV(~  y2  -  kj) |  **  1  in 
the  range  of  values  of  x  for  which  ifi(x)  and  if,  are 
appreciably  different  from  zero.  This  allows  us  to  expand 
the  radical  in  (37)  using  the  binomial  theorem  and  retain¬ 
ing  only  a  finite  number  of  terms.  This  results  in 


/«'(y)  =  (-V-*i),/2j/  jh'fidx 


2(~Y  2~kl) 


/» _ 

XMrdX 

- 


_ 1 

8  (-T*-*.’)1 

•/_*X4M,«/X+  •••)•  (38) 

The  integrals  do  not  depend  on  y,  so  they  need  to  be 
computed  only  once  for  a  given  set  of  material  parameters, 
and  not  for  every  value  of  y,  as  (31)  would  require. 
Furthermore,  the  properties  of  the  Fourier  transform 
guarantee  that 

/“  MrdX  *=«//••  (39) 

•'-00 

Agreement  with  exact  computation  of  (31)  is  excellent  with 
the  first  three  terms  in  the  ^-pansion  (38)  and  the  compu¬ 
tation  time  is  reduced  snhstnn*ially.  by  nearly  two  orders 
of  magnitude.  This  allows  us  to  perform  more  extensive 
modeling  using  this  method. 

Fig.  7  shows  plots  of  the  modal  loss=Re(y)  as  a 
function  of  An,  with  g0,  the  gain  under  the  stripe,  as  a 
Parameter.  We  notice  that  for  each  value  of  g0,  An  can  be 
decreased  up  to  a  certain  value  beyond  which  the  funda¬ 
mental  mode  becomes  lcatcy,  i.e.,  the  antiguiding  effect  of 
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Fig.  6.  Mode  attenuation  Re(y)  and  number  of  discrete  trapped  modes 
versus  An. 


the  (negative)  An  offsets  the  guiding  effect  of  the  gain 
distribution. 

More  negative  values  of  A  n  are  required  to  offset  higher 
values  of  g0.  For  any  given  g0,  using  (5),  (18),  and  (20)  with 
/  =  0,  it  can  be  shown  that  the  boundary  value  of  An  is 

An  =  -ix5n0(g0  +  as)2.  (40) 

These  values  are  marked  with  vertical  lines  in  Fig.  7.  This 
relation  shows  that  the  guiding  effect  of  the  gain  distri¬ 
bution  does  not  depend  on  the  individual  values  of  g0  (gain 
under  the  stripe)  and  -  as  (gain  far  away  from  the  stripe), 
but  only  on  their  difference  g0  —  (  -  as)  =  g0  +  as. 

V.  Effective-Index  Calculation 

The  effective-index  method  consists  basically  in  reducing 
a  two-dimensional  problem  to  an  equivalent  one-dimen- 
sional  one.  In  our  case,  the  two-dimensional  character  of 
the  problem  is  given  by  the  dependence  of  the  dielectric 
constant  on  x  and  y.  As  a  first  approximation,  the  varia¬ 
tion  in  one  direction  (in  our  case:  x)  is  neglected;  this  is 
justified  if  this  variation  is  much  less  than  that  in  the  y 
direction.  This  is  equivalent  to  approximating  the  wave¬ 
guide  with  a  simple  three-layer  guide  whose  dielectric 
constants  do  not  vary  with  x.  The  solution  of  this  problem 
yields  the  transverse  or  vertical  variation  of  the  field.  Next, 
the  original  equation  describing  the  two-dimensional  equa¬ 
tion  in  x  can  be  solved  for  the  lateral  variation  of  the  field, 
and  the  overall  solution  is  approximated  by  the  product  of 
this  lateral  field  and  the  vertical  field  found  from  the 
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.  three-layer  problem.  We  start  with  the  wave  equation  in 
two  dimensions,  which  is  obtained  merging  (7)  and  (8) 

V,j*  +  [y2  +  *oM*.>)]*  =  0  (41) 

where 


*(x,y)  = 


'Ay 


tr  + 


Ak 


cosh2(x/x0)  ’ 


\y\>d/2 
l>’l  <  d/2 ' 


For  the  simple  three-layer  guide  we  assume  k(*): 
inside  the  active  layer.  We  then  have 


d2<(> 
dy 2 


P2<t> 


with 


p2  +  q2  =  kl[Ko-icA] 
/>  =  <)  tan(<j<//2). 


'-d/2 


/  <t>*{y)4>(y)dy 

»  —  OVN 


We  obtain 

f  . 

d24  , 

i 

dx1 

y2  -  qlt  +  *o*s,ff  +  *o 


coshJ(x/x0) 


with 


?«ff  =  r*o*o  ~  P2  ~  *0*4 


KSt(l  -  ^*s 

Akc((  =  FA*c. 


\y\*kd/2  (43a) 

ly\>d/2  (43b) 

(44a) 
(44b) 

Now,  we  transform  (41)  making  ^(x,  y)  =  ^(•x)<f>(>').  mul¬ 
tiplying  it  by  <#>*(>•),  and  integrating  it  overy  from  -  oo  to 
oo,  taking  into  account  (42)  and  (43)  and  defining  the 
confinement  factor 

jd/ 1  't>*(y)4>(y)dy 


Fig  8. 


Mode  attenuation  Re{y)  versus  An  for  numeric  and  effective-  •> 
index  methods,  for  g0  -  200  cm  “ 1 .  ; 


(45) 

=  0 

(46) 

(47a) 
(47b) 
(47c) 

Equation  (46)  has  the  same  form  as  (10b).  It  will  also  have 
polynomial  solutions  similar  to  (17)  that  represent  trapped 
modes 

+,(x)  =  [cosh  ( x/x0 )] '  ‘  * /+  ‘I  ( tanh  j- ) 

/  =  0, 1,2,  •  •  •  (48) 

where  60tff,  Aicef(,  and  kScU  satisfy  relations  identical  to 
(12b).  (15b),  and  (18).  For  the  fundamental  mode,  /  =  0, 
B  =  h0t(f,  and  using  (44b),  we  obtain 

72“-*0*4-%I-(p2-*0rAK)  (49) 

*o 

for  the  propagation  constant. 

VI.  Discussion 

Values  of  y  obtained  using  (49)  are  compared  with  those 
obtained  with  the  numerical  method  in  Fig.  8.  Solutions 
are  very  close  for  all  values  of  An  for  which  the  mode 


I 

!» 


Fig.  9.  Peak  active  region  gain  g0  versus  An  for  severat  values  of  Re(y\ 
Dotted  area  indicates  region  where  fundamental  mode  becomes  leaky 


exhibits  a  gain  which  is  relatively  high  and  with  !o  - 
sensitivity  to  An.  The  results  differ  most  in -the  range  of .'  i 
for  which  the  mode  has  a  net  loss  or  has  a  relatively  hv. 
gain  with  higher  sensitivity  to  An.  This  is  the  same  region 
for  which  only  one  trapped  mode  exists,  so  the  error  in  »V 
numerical  method  is  the  greatest  because  the  neglected 
continuum  is  more  important. 

Fig.  9  gives  the  required  value  of  the  peak  power  gain  g0 
under  the  stripe  to  obtain  a  given  modal  gain  C,  as  :> 
function  of  An,  using  the  effective-index  method.  Also 
included  is  the  region  for  which  the  fundamental  mode 
becomes  leaky.  Figs.  10-12  show  the  normalized  lateral 
field  distributions  for  different  values  of  An.  The  increas¬ 
ing  antiguiding  effect  of  decreasing  An  is  apparent.  The 
distance  x  at  which  the  gain  is  zero  (loss/gain  boundary)  is 
shown  in  these  figures  with  vertical  dashed  lines.  It  can  be 
shown  to  be  given  by 

"(,+5-?r-  <5o) 


x  =  x0cosh" 


For  the  set  of  parameters  considered.  3f=l.45  x0.  This 
allows  us  to  qualitatively  understand  why  the  modes  have 
the  loss  (gain)  indicated.  The  vertical  confinement  factor  T 
did  not  vary  appreciably  with  An;  a  typical  value  for  the 
case  considered  was  T  ~  0.4963.  We  see  that  lateral  varia¬ 
tions  in  the  refractive  index  affect  the  gain  much  more  by 
altering  the  lateral  field  distribution  than  by  affecting  the 
vertical  variation. 


/f  is,  *’■  .  ”  - 


VT 


•X 

l 


A 


c 


Fig.  10.  Normalized  lateral  field  distribution,  for  An-  —  0.01 56.  The 
resulting  modal  loss  is  20  cm'1.  Vertical  dashed  lines  indicate  the 
distance  at  which  the  gain  inside  the  active  layer  is  zero. 


Fig.  II.  Normalized  lateral  field  distribution  for  An  = -0.006.  Re(y)- 
-4.46  cm" 1  (gain).  Vertical  dashed  lines  have  the  same  meaning  as  in 
Fig  10. 


LATERAL  DISPLACEMENT  X/Xo 

Fig  12.  Normalized  lateral  field  distribution  for  An  —  0.002.  The  mods! 
gain,  -  Ref  y),  is  now  32.3  cm' 

j  The  effective-index  method  is  seen  to  be  a  fast  and 
!  relatively  accurate  way  to  obtain  the  field  distributions  for 
|  the  class  of  waveguides  considered,  for  which  the  numerical 
method  we  used  is  not  practical  for  extensive  modeling  due 
to  its  long  computation  time  in  spite  of  all  approximations. 
The  remarkable  agreement  between  the  effective-index 
I  method  and  experimental  results  found  by  other  workers 
.  (|  101)  increases  our  confidence  in  this  powerful  approxi- 
:  mate  method. 
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APPLICATION  OF  A  GENERAL  PURPOSE  CIRCUIT  SIMULATION  PROGRAM  TO 
SEMICONDUCTOR  LASER  MODELING  -  PART  I  s  ELECTRICAL  ANALYSIS* 

A.  Linz  and  J.K.  Butler 
Electrical  Engineering  Department 
Southern  Methodist  University 
Dallas,  Texas  75275 

Abstract 

This  article  describes  the  application  of  SPICE  to  the  electrical 
modeling  of  modern  semiconductor  laser  structures  which  cannot  be  described 
analytically.  Given  the  geometrical  structure  of  the  device  and  some  material 
characteristics  such  as  impurity  concentration  and  type  of  semiconductor  (p  or 
n),  we  must  obtain  the  distribution  of  potential  and  current  density.  We  are 
mainly  interested  in  knowing  the  current  density  at  the  edge  of  the  active 
layer,  since  it  determines  the  distribution  of  injected  carriers  inside  it, 
together  with  the  diffusion  equation.  The  carrier  concentration,  in  turn, 
determines  the  optical  gain  in  the  active  layer.  The  gain  profile  will  have  a 
strong  influence  on  the  lateral  dependence  of  the  optical  field.  In  this 
article  we  start  analyzing  a  very  simple  structure  which  nevertheless  contains 
many  features  found  in  real  lasers.  Then  we  present  a  discrete  model  which  can 
be  applied  to  virtually  any  geometry  and  we  use  it  to  explore  some  laser 
structures  with  nonlinear  elements  (p-n  junctions).  We  confine  ourselves  to 
purely  electrical  effects;  the  interaction  between  carrier  diffusion  and 
electrical  parameters  is  left  for  a  subsequent  article. 

*  Supported  in  part  by  the  U.  S.  Army  Research  Office  and  RCA  Laboratories. 


APPLICATION  OF  A  GENERAL  PURPOSE  CIRCUIT  SIMULATION  PROGRAM  TO 
SEMICONDUCTOR  LASER  MODELING  -  PART  II  :  DIFFUSION  ANALYSIS* 

A.  Linz  and  J.K.  Butler 
Electrical  Engineering  Department 
Southern  Methodist  University 
Dallas,  Texas  75275 

Abstract 

This  article  is  the  continuation  of  a  previous  one  in  which  the  modeling 
of  electrical  characteristics  of  semiconductor  laser  devices  using  SPICE  was 
described.  Here  we  consider  the  application  of  SPICE  to  the  modeling  of 
diffusion  and  the  integration  of  the  electrical  and  the  diffusion  models.  The 
model  is  then  applied  to  a  particular  structure  and  the  results  are  discussed. 
The  model  shows  that  the  junction  current  density  distribution  remains 
virtually  unchanged  for  a  wide  range  of  diffusivities. 

*  Supported  in  part  by  the  U.  S.  Army  and  RCA  Laboratories. 


AN  INTEGRATED  ELECTRICAL/DIFFUSION  MODEL 
FOR  THE  P-N  JUNCTION* 

A. Linz,  J.K.  Butler,  and  K.H.  Heizer 
Electrical  Engineering  Department 
Southern  Methodist  University 
Dallas,  Texas  75275 

Abstract 

A  lumped  parameter  circuit  model  is  presented,  which  solves  for  the 
coupled  electrical/diffusion  problem  for  a  p-n  junction.  Externally,  the 
model  yields  the  electron  and  hole  current  components  separately.  Internally, 
an  adjoint  "diffusion  network"  which  is  an  electrical  equivalent  of  the  one¬ 
dimensional  diffusion  equation,  yields  the  values  of  excess  carrier 
concentration  in  both  the  n  and  p  regions.  The  model  is  then  used  to  obtain 
the  transient  response  of  a  diode.  The  effect  of  the  material  parameters 
(lifetime,  diffusivity,  resistivity)  on  the  electrical  response  of  the  device 
can  be  readily  modeled.  When  used  with  SPICE,  the  model  can  easily  be  used  as 
a  subcircuit  to  incorporate  it  into  networks  containing  diodes. 

*  Supported  in  part  by  the  U.  S.  Army  Research  Office  and  RCA  Laboratories. 
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The  lateral  dielectric  constant  in  stripe  geometry  injection 
lasers  is  modeled  with  a  step  and  a  parabolic  dependence.  For  the 
step  profile  a  set  of  cutoff  curves  for  the  lateral  modes  is  used 
to  compare  various  laser  designs.  As  the  vertical  geometry  of  a 
laser  is  varied  a  contour  is  described  on  the  cutoff  curves  which 
is  a  measure  of  lateral  mode  stability.  Two  laser  designs  examined 
are  the  simple  stripe  geometry  and  channeled-substrate  planar  lasers. 
The  parabolic  model  of  the  lateral  dielectric  constant  is  extended 
to  include  matched  fields  at  the  vertical  heterojunctions.  A  result 
is  a  deformed  fundamental  mode  shape  in  the  case  of  small  normalized 
frequencies  and  depression  of  the  refractive  index  under  the  stripe. 
This  presents  an  alternative  explanation  for  the  shoulders  that  have 
heretofore  been  attributed  to  leaky  waves. 

Lateral  mode  stability  with  drive  current  above  threshold  is 
complicated  by  an  uncertain  dielectric  profile  perturbed  by  free 
carrier  effects  and  spatial  hole  burning.  An  experimental  scheme 
is  devised  to  measure  the  free  carrier  effect  on  the  refractive 


index  in  GaAs  active  regions  of  injection  lasers.  Transverse  field 
and  modal  selection  mechanisms  of  the  active  cavity  are  character¬ 
ized  for  a  broad  area  high  power  laser  as  a  function  of  temperature 
and  injection  level.  The  amount  of  index  depression  as  a  result  of 
carrier  increase  in  the  active  region  is  measured  to  be  (3  -  9)  x 
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A  set  of  programs  to  model  the  operation  of  semicon¬ 
ductor  laser , devices  below  and  above  threshold  Is  presen¬ 
ted,  together  with  the  design  and  construction  of  an  auto¬ 
matic  measurement  system  used  to  obtain  far  field  radiation 
patterns.  The  programs  incorporate  electrical  ,  diffusion 
and  optical  modeling. 

The  program  SPICE  is  used  to  model  electrical  and 
diffusion  phenomena,  while  direct  numerical  techniques  are 
used  to  solve  the  optical  and  diffusion  problems  . 
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